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Introduction 

Many  physical  mechanisms  can  drive  a  laser  from  a  single-mode  to  a  multi-mode 
regime.  Common  examples  are  spatial  and  spectral  hole  burning  (SHB),  saturable 
absorption,  and  self-phase  modulation  [1,2,3],  Understanding  these  mechanisms  is  of 
key  importance  to  laser  science  and  technology,  whether  one  is  interested  in 
single-mode  behavior  or  in  a  particular  multi-mode  operation  such  as  mode  locking. 
The  multimode  regimes  listed  above  are  well  understood  and  documented  both 
theoretically  and  experimentally,  in  many  types  of  lasers.  However  the  understanding 
of  multimode  regimes  in  quantum  cascade  lasers  (QCLs)  is  still  in  its  infancy,  as  these 
lasers  were  only  demonstrated  in  1994  [4],  and  studies  of  their  multimode  regimes 
commenced  even  more  recently  [5, 6, 7, 8], 

As  it  was  shown  recently  [9]  that  multimode  dynamics  in  QCLs  are  different  from  that 
of  more  common  lasers.  This  is  mainly  due  to  the  unusually  fast  gain  recovery  of  QCLs, 
which  occurs  on  a  picosecond  scale.  While  a  saturable  absorber  triggers  mode  locking 
in  lasers  with  slow  gain  recovery  (relative  to  the  roundtrip  time),  in  lasers  with  fast  gain 
recovery  a  saturable  absorber  triggers  a  mechanism  similar  to  the 
Risken-Nummedal-Graham-Haken  (RNGH)  instability  [10,11].  While  in  standard 
semiconductor  lasers  carrier  diffusion  eliminates  spatial  hole  burning,  in  QCLs  the  gain 
recovery  process  is  faster  than  carrier  diffusion,  and  spatial  hole  burning  is  dominant. 
In  this  report  we  present  a  detailed  study  of  multimode  regimes  in  QCLs.  The  results  of 
Ref.  [9]  are  substantiated  and  extended.  The  first  part  of  the  paper  is  theoretical  and  the 
second  one  is  experimental.  In  the  theoretical  section,  the  laser  Maxwell-Bloch 
equations  are  introduced  and  analyzed  for  a  Fabry-Perot  cavity.  This  way  one  can 
study  the  interplay  of  coherent  phenomena  and  spatial  hole  burning.  A  saturable 
absorber  is  added  to  the  model  as  well.  The  model  is  studied  analytically  and 
numerically.  The  stability  region  of  a  continuous  wave  (CW)  solution  is  found.  It  is 
shown  that  in  a  ring  cavity,  the  presence  of  a  saturable  absorber  lowers  the  threshold  of 
the  RNGH  instability  from  about  nine  times  above  laser  threshold  to  arbitrarily  low 
above  laser  threshold,  depending  on  the  strength  of  the  absorber.  However,  the  nature 
of  the  instability  remains  the  same:  the  population  inversion  begins  to  oscillate  at  the 
Rabi  frequency,  modulating  the  gain  in  the  laser.  The  result  is  sidebands  around  the 
original  CW  mode,  separated  from  it  by  roughly  the  Rabi  frequency.  A  Rabi  splitting  in 
the  spectrum  is  the  primary  signature  of  the  RNGH  instability. 

It  is  then  shown  that  the  ring-cavity  Maxwell-Bloch  model  with  a  saturable  absorber 
can  explain  the  Rabi  splitting  and  the  lowering  of  the  threshold,  but  cannot  explain  a 


key  feature  in  the  experimental  spectra.  The  latter  exhibit  two  relatively  equal  groups  of 
modes  with  a  gap  in  between,  whereas  the  ring-cavity  Maxwell-Bloch  model  predicts  a 
large  central  mode  in  the  spectrum  with  two  sidebands.  In  view  of  this  discrepancy,  the 
Maxwell-Bloch  equations  were  extended  to  include  coupling  between 
counter-propagating  modes  in  a  Fabry-Perot  cavity,  which  is  of  course  a  more  adequate 
model  for  QCLs.  This  model  takes  into  account  the  development  of  SHB.  The  inclusion 
of  SHB  generates  theoretical  spectra  with  a  Rabi  splitting  and  without  a  central  peak  in 
the  spectrum,  in  agreement  with  experiments. 

After  presenting  the  theoretical  results,  a  thorough  study  of  the  experimental 
phenomenology  is  presented.  It  is  shown  that  QCLs  with  narrower  active  regions  tend 
to  exhibit  a  more  pronounced  Rabi  splitting  than  lasers  with  wider  active  regions. 
Lasers  with  a  wider  active  region  tend  to  exhibit  multimode  spectra  that  are  governed 
by  SHB.  The  explanation  we  propose  to  this  behavior  is  that  narrow  QCLs  have  a 
stronger  saturable  absorption  effect  than  wider  ones  as  a  result  of  a  Kerr-lensing  effect 
[3,5],  In  this  case  a  nonlinear  index  enhancement  in  the  waveguide  core  gives  rise  to  an 
increased  overlap  of  the  transverse  laser  mode  with  an  active  region  and  a  reduced 
overlap  with  lossy  sidewalls,  leading  to  an  additional  enhancement  of  the  saturable 
absorption.  The  idea  of  saturable  absorption  by  Kerr-lensing  is  further  supported  by  the 
fact  that  QCLs  with  metal  coating  on  the  sides  of  the  ridge  have  a  stronger  RNGH 
behavior  than  lasers  without  metal  coating:  The  metal  coating  enhances  saturable  losses 
originated  from  Kerr-lensing. 

In  the  last  part  of  the  paper,  the  temperature  dependence  of  the  multimode  behavior  in 
QCLs  is  studied  experimentally.  It  is  found  that  at  higher  temperature  the  multimode 
behavior  tends  to  be  governed  by  the  RNGH  instability,  whereas  at  lower  temperatures 
it  is  governed  by  SHB.  This  behavior  can  be  explained  by  the  fact  that  at  higher 
temperatures  hot  carriers  populate  more  states  in  the  injector  superlattice,  creating 
aditional  quasi-resonant  absorption  transitions  between  ground  and  excited  minibands. 
This  leads  to  an  additional  saturable  absorption  of  laser  radiation. 

This  paper  is  organized  as  follows:  Section  2  gives  a  brief  survey  of  prior  work  on 
multimode  regimes  in  QCLs,  Sections  3  study  theoretically  the  Maxwell-Bloch 
equations  in  a  Fabry-Perot  cavity,  Sections  4-6  summarize  the  experimental  study,  and 
Section  7  is  a  brief  summary. 

2.  Prior  work  on  multimode  regimes  in  QCLs 

Multimode  regimes  in  QCLs  were  observed  in  a  series  of  recent  works  [5, 6, 7, 8],  In 
Ref.  [5],  for  example,  it  was  observed  that  at  a  certain  pumping  current  above  lasing 
threshold,  QCLs  cease  to  operate  in  CW  and  develop  a  multimode  regime.  This 
multimode  regime  was  characterized  by  a  broadband  optical  spectrum  and  a  narrow 
(less  than  1 00  kHz)  radio  frequency  (RF)  beat  note  in  the  power  spectrum. 

The  narrow  beatnote,  whose  width  is  1/105  of  the  central  frequency,  shows  that  the 
waveform  of  the  electric  field  circulating  in  the  laser  cavity  was  stable  over 
approximately  105  roundtrips.  In  other  words,  the  phase  relationships  between  the 
longitudinal  modes  were  stable  for  about  105  roundtrips.  The  modes  were  therefore 
locked.  However  in  order  to  characterize  the  waveform  which  is  circulating  in  the  laser 
and  to  see  if  it  is  indeed  an  isolated  pulse,  as  in  traditional  mode  locking,  one  has  to 
apply  pulse  characterization  techniques,  such  as  second  order  autocorrelation. 

At  the  time  when  the  experiment  in  Ref.  [5]  was  performed,  no  second-order 
autocorrelation  apparatus  was  available.  However  second  harmonic  generation  from 
QCLs  [6]  provide  some  information  for  pulse  characterization.  This  measurement 


shows  an  increase  by  more  than  a  factor  of  5  in  the  second  harmonic  signal  as  the 
multimode  behavior  sets  in.  This  increase  indicates  that  the  duty  cycle  of  the  pulses  was 
roughly  5.  Since  the  number  of  modes  in  the  spectrum  is  more  than  5,  we  conclude  that 
not  all  the  modes  were  locked  into  a  pulse.  Without  better  pulse  characterization  data, 
one  can  infer  that  the  laser  could  have  had  more  than  one  pulse  per  roundtrip. 

It  should  also  be  noted  that  traditional  mode  locking,  with  a  single  pulse  per  roundtrip, 
requires  that  the  gain  recovery  time  be  longer  than  the  cavity  roundtrip  [3],  In  QCLs 
this  condition  is  violated,  and  therefore  according  to  the  existing  theory  one  cannot 
expect  mode  locking  with  a  single  pulse  per  roundtrip  in  QCLs.  One  can  expect 
multiple  pulses  per  roundtrip. 

In  view  of  the  above,  one  can  see  that  the  nature  of  various  multimode  regimes  in  QCLs 
requires  further  elucidation.  This  is  the  main  goal  of  the  present  work. 

3.  Maxwell-Bloch  equations  in  a  Fabry-Perot  cavity 
Theory 


Numerical  studies  based  on  the  master  equation  of  mode  locking  [2,  3]  with  a  1.5  ps 
gain  recovery  time  (which  represents  the  relaxation  time  of  electrons  between  the  states 
of  the  laser  transition)  did  not  lead  to  stable  self  modelocking.  Physically  for  the  latter 
to  occur  the  gain  recovery  time  must  be  greater  than  the  roundtrip  time  to  prevent  the 
formation  of  pulses  at  intervals  smaller  than  the  latter  The  mechanism  responsible  for 
the  partial  mode  locking  observed  in  the  experiments  was  therefore  yet  to  be 
understood. 

During  the  past  year,  we  have  developed  various  phenomenological  models  that  might 
explain  the  currently  observed  multimode  dynamics  and  hopefully  give  guidance  on 
how  to  achieve  complete  mode  locking  of  QCLs.  The  most  consistent  picture  comes 
from  a  model  which  includes  a  saturable  absorber  and  which  is  based  on  the 
Maxwell-Bloch  equations  for  a  two-level  active  medium.  Although  the  QCL  gain 
medium  should  be  described  by  a  more  complicated  set  of  equations,  the  main  physics 
can  already  be  understood  from  this  theory. 

The  Maxwell-Bloch  equations  are  well  known  to  exhibit  a  modulation  instability  first 
reported  by  Risken  and  Nummedal  and  independently  Graham  and  Haken  [4,  5].  The 
mechanism  responsible  for  the  onset  of  the  Risken-Nummedal-Graham-Haken  (RNGH) 
instability  is  related  to  coherent  Rabi  oscillations  of  population  inversion  on  the  laser 
transition  that  give  rise  to  sidebands  separated  by  the  Rabi  frequency  from  the  initial 
single-frequency  cw  emission.  The  instability  develops  only  when  the  laser  power  is 
high  enough  so  that  the  Rabi  frequency  exceeds  the  relaxation  rate  of  the  population 
inversion.  For  lower  frequencies  the  Rabi  oscillations  of  inversion  are  overdamped  and 
do  not  lead  to  the  instability.  The  appearance  of  sidebands  in  the  spectrum  is  manifested 
by  high-frequency  pulsations  in  the  time  domain.  The  dynamics  far  above  the 
instability  threshold  can  be  very  complicated  and  rich  -  from  regular  mode-locking  to 
completely  chaotic  pulsations. 

This  mechanism  can  explain  the  main  features  of  the  measurements  described  in  the 
experimental  part.  The  RNGH  instability  requires  in  principle  very  high  pumping 
powers,  typically  9  to  10  times  above  threshold.  In  our  model  however,  the  threshold 


for  the  instability  can  be  lower  by  more  than  a  factor  of  four,  because  of  the  presence  of 
the  saturable  absorber.  This  result  can  be  established  analytically  by  linear  stability 
analysis,  as  well  as  by  numerical  simulations.  The  threshold  for  the  instability  is  thus 
easily  reachable  in  the  experiment.  Saturable  absorption  in  our  devices  is  due  to  the 
presence  of  metal  on  the  side  of  the  laser  ridge  (see  Fig.  1  (b)).  It  is  the  source  of 
non-negligible  waveguide  losses,  which  decrease  at  large  intensities  thanks  to  Kerr 
lensing  [1],  Our  future  plans  include  the  design  of  a  intra-cavity  saturable  absorber  due 
to  resonant  intersubband  absorption,  which  may  be  easier  to  control  than  the  losses  in 
sidewalls. 

Note  that  since  short  time  scales  comparable  to  the  ultrashort  gain  recovery  time  are 
important,  taking  into  account  the  full  phase  relationship  between  the  medium 
polarization  and  the  optical  field  is  crucial.  Indeed  the  master  equation  of  mode  locking, 
in  which  the  polarization  is  adiabatically  eliminated,  gave  no  pulses  at  all  for  the  same 
set  of  parameters. 

Figure  5  (a)  shows  that  above  a  certain  threshold  power,  trains  of  pulses  as  short  as  the 
gain  recovery  time  (1.5  ps)  are  obtained  in  our  simulations.  The  periodicity  in  the  time 
domain  translates  into  a  corresponding  modulation  that  can  be  seen  in  the  spectrum 
shown  in  Fig.  5  (b).  This  periodicity  can  be  identified  with  oscillations  of  the 
population  inversion  at  the  Rabi  frequency.  As  the  laser  field  becomes  stronger,  the 
separation  between  the  different  groups  of  modes  in  the  optical  spectra  grows. 

The  parameters  used  in  the  simulations  displayed  in  Fig.  5  correspond  to  the  laser 
characteristics  and  the  experimental  conditions  of  Fig.  2  (a).  All  the  parameters  were 
measured  independently,  apart  from  the  saturable  absorber  coefficient.  The  latter  was 
calculated  from  the  measured  instability  threshold.  The  spectral  separation  between  the 
two  groups  of  modes  in  Fig.  2  (a)  is  fairly  close  (within  a  factor  of  1.5)  to  the  separation 
in  the  simulated  spectra.  Although  a  few  parameters  such  as  the  saturation  intensity  of 
the  saturable  absorber  can  only  be  guessed,  our  model  is  consistent  with  our 
experimental  data  to  a  very  good  extent. 

Another  evidence  comes  from  the  direct  comparison  between  the  experimental  splitting 
and  the  Rabi  frequency,  which  can  be  calculated  from  the  measured  field  intensity  and 
the  dipole  matrix  element.  For  an  output  power  of  lOOmW,  corresponding  to  the  trace 
at  a  current  of  700  mA  in  Fig.  2  (a)  the  Rabi  frequency  is  0.7  THz,  which  is  in  excellent 
agreement  with  the  experimental  separation  between  the  group  of  two  modes  (0.8THz). 
The  corresponding  Rabi  angular  frequency  (2n  times  the  Rabi  frequency)  is 
considerably  greater  than  the  inverse  of  the  gain  recovery  time,  thus  satisfying  the 
instability  criterion. 

There  is  however  an  important  qualitative  difference  between  experiment  and  theory: 
the  latter  always  has  a  large  central  peak  in  the  spectrum  with  sidebands  (groups  of 
modes)  of  similar  intensity,  while  the  former  is  not  visible  in  the  experiment 
consistently  shows  only  two  groups  of  modes  with  similar  intensity.  Our  numerical 
studies  of  the  Maxwell-Bloch  equations  shows  that  such  solutions  do  not  occur.  It  is 
therefore  yet  to  be  understood  which  additional  physical  mechanism  should  be  included 
in  the  theoretical  model  in  order  to  eliminate  this  discrepancy. 


Note  that  the  emission  of  the  pulses  becomes  increasingly  chaotic  as  a  function  of  the 
pumping  power.  Our  simulations  also  predict  a  modulation  of  the  laser  intensity  at  the 
cavity  round-trip  time  as  observed  experimentally  with  a  fast  QWIP  and  a  spectrum 
analyzer.  This  modulation  shows  that  there  is  a  non-trivial  phase  relation  between  the 
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Figure  5:  Results  of  the  simulations,  gain  recovery  time  1.5  ps.  (a)  shows  the  intensity  as 
function  of  time  during  10  roundtrips,  and  (b)  shows  the  optical  spectrum.  The  number 
on  the  left  of  each  subplot  is  the  pumping  current  divided  by  the  lasing  threshold  current. 
The  intensity  of  the  peak  at  zero  frequency  is  normalized  to  1  because  it  is  much  stronger 
than  the  sidebands. 


longitudinal  modes,  or  equivalently  between  the  pulses.  This  aspect  of  the  simulations 
is  is  important  as  it  may  help  to  understand  the  shape  of  the  autocorrelation  trace. 


4.  Experimental  results:  ridge  lasers 

We  first  study  the  multimode  regimes  in  standard  ridge  QCLs,  in  which  the  sidewalls  of 
the  laser  ridges  are  covered  by  a  thick  layer  of  electrically  plated  gold  contact.  This  acts 
as  a  Kerr-lens  type  saturable  absorber;  see  Fig.  8  for  a  typical  laser  cross  section.  The 


active  region  of  the  samples  tested  is  based  on  a  three-quantum-well  design  emitting  at 
a  wavelength  X~8pm  [31].  The  wafer  was  grown  by  metalorganic  vapor  phase  epitaxy 
(MOVPE).  Fig.  9  a  shows  the  voltage-current  (V-I)  and  light-current  (L-I) 
characteristics  of  a  10  pm  wide  laser  operated  in  CW  at  200K,  and  b  shows  the 
corresponding  optical  spectra.  The  laser  was  cleaved  into  a  2  mm  long  bar  and  soldered 
with  Indium  onto  a  copper  heat  sink.  The  optical  power  was  measured  by  an  OPHIR 
thermal  head  powermeter  with  a  collection  efficiency  of  nearly  100%.  The  spectra  were 
measured  by  a  Nicolet  Fourier  transform  infrared  spectrometer  (FTIR)  equipped  with  a 
deuterated  triglycine  sulphate  (DTGS)  detector. 


8.  Scanning  electron  microscope  (SEM)  image  of  the  cross  section  of  a  10  [Trial  modejm  wide 

ridge  QC  laser. 

o 


FIG.  9:  (a)  V-I  &  L-I  curves  and  (b)  optical  spectra  vs.  pumping  ratio  (j/jtt)  above  threshold  obtained  in  CW  at  200K  of 
a  10  pm  wide  ridge  laser  emitting  at  8  pm  (wafer  #  2721)  (c)  Spectral  splitting  and  twice  the  Rabi  frequency  fiRabi/(2ir)  vs. 
square  root  of  output  power  collected  from  a  single  laser  facet.  The  different  quantities  reported  on  the  graph  were  deduced 
from  the  experimental  data  shown  in  9  (b).  The  dashed  line  is  a  least-square  linear  fit  of  the  data. 


As  shown  in  Fig.  9  b,  the  laser  spectrum  is  single  mode  close  to  laser  threshold.  It 
broadens  and  splits  into  two  separated  humps  as  the  pumping  current  increases.  The 
separation  between  the  two  peaks  of  the  two  humps  increases  linearly  with  the  square 


root  of  the  collected  output  power  from  one  facet,  as  shown  in  Fig.  9  c.  The  Rabi 
angular  frequency  can  be  calculated  from  the  collected  output  power,  using  the  formula 
Q  Rahi  =  fiE/ti,  where  p  is  the  electron  charge  times  the  matrix  element  of  the  laser 

transition  (=1.9  nm  for  this  particular  device).  Iave  is  the  average  intracavity  intensity  in 
the  gain  region,  which  can  be  derived  from  the  measured  output  power  [32].  For  all 
values  of  intensity  corresponding  to  the  spectra  reported  in  Fig.  9  b,  D.Rahi/2 n  was 
calculated,  multiplied  by  a  factor  of  two  and  then  added  to  Fig.  9  c  (solid  line). 
Reasonably  good  agreement  is  found  between  the  experimental  splitting  and  twice  the 
estimated  Rabi  frequency.  The  error  bars  of  the  spectrum  splittings  come  from  the 
uncertainty  in  determining  the  exact  position  of  the  peaks,  which  is  the 
full-width-at-half-maximum  (FWHM)  of  the  humps.  As  mentioned  previously  in  the 
theoretical  section,  the  RNGF1  instability  predicts  that  large  intracavity  intensity  will 
result  in  parametric  gain  at  frequencies  detuned  from  the  maximum  of  the  gain  curve  by 
the  Rabi  frequency.  The  measured  spectra  thus  show  strong  indication  of  the  RNGH 
instability  in  ridge  QCLs. 
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FIG.  10:  (a)  The  V-I  V  L-I  curve  and  (b)  the  optical  spectra  vs.  pumping  ratio  (j/jih)  above  threshold  obtained  from  a  10 
ni  inn  wide  ridge  QCL  operated  in  CW  at  77K  emitting  at  8/»m  (wafer  #  2743).  Note  that  mode  hopping  occurs  while  there  is 
a  jump  in  the  L-I  curve. 


The  lowering  of  the  RNGF1  instability  threshold  in  our  QCLs  is  due  to  the  presence  of  a 
saturable  absorber.  This  phenomenon  is  demonstrated  analytically  in  the  theoretical 
section.  Such  a  saturable  absorption  mechanism  in  our  experiments  is  likely  to  come 
from  Kerr-lensing,  caused  by  a  nonlinear  (i.e.  intensity  dependent)  refractive  index  ml 
in  the  active  region  [5],  As  the  light  intensity  increases,  the  mode  becomes  more 
confined  in  the  plane  transverse  to  the  propagation  direction,  and  the  net  gain  it 
undergoes  also  increases.  The  reason  is  twofold:  First,  the  mode  overlaps  more  with  the 
active  region,  leading  to  a  larger  modal  gain  (this  mechanism  is  often  called  "soft 
Kerr-lensing").  Second,  the  overlap  with  the  metal  contacts  is  reduced,  leading  to 
smaller  losses. 


11.  The  experimental  setup  of  a  two-photon  autocorrelation  measurement.  (Inset)  Conduction 
band  diagram  of  the  two-photon  QW1P  showing  three  equidistant  energy  levels. 

The  same  RNGH  splitting  in  spectra  is  observed  in  many  different  devices,  from  wafers 
grown  by  both  molecular  beam  epitaxy  (MBE)  and  MOVPE.  Fig.  10  shows  the  spectra 
of  a  laser  fabricated  from  a  MBE-grown  wafer  with  the  same  active  region  design  as  in 
Fig.  7,  taken  in  continuous  wave  at  77K.  The  spectrum  starts  from  single  mode  close  to 
the  laser  threshold,  and  the  mode  hops  at  a  pumping  ratio  of  j/jth  =1.66,  where  there  is  a 
correponding  kink  in  the  L-I  curve  at  330mA.  At  590mA,  there  is  another  kink  in  the 
L-I  curve,  relating  to  the  broadening  of  the  spectrum  and  a  hop  in  the  mode  center. 
Mode-hopping  and  its  relationship  to  kinks  in  the  L-I  curves  are  common  phenomena  in 
semiconductor  lasers  [33,34,35],  but  its  cause  in  QC  lasers  is  not  yet  studied  and  is 
beyond  the  focus  of  this  paper.  After  the  spectrum  broadens,  it  forms  two  separated 
humps  whose  peaks  shift  apart  with  increased  optical  power,  similar  to  the 
MOVPE-grown  case. 

Characterization  of  our  ultra-short  pulses  was  done  using  the  standard  method  of 
second-harmonic  interferometric  autocorrelation.  The  setup  is  based  on  a  Michelson 
interferometer  in  which  the  input  beam  is  split  into  two  and  one  of  them  is  delayed  by 
r.Once  recombined,  the  two  pulses  are  sent  colinearly  first  into  a  nonlinear  crystal,  and 
then  a  filter,  which  allows  only  the  second-harmonic  generation  (SHG)  component  to 
be  detected.  One  can  test  if  it  is  an  isolated  pulse  from  the  ratio  between  the  interference 
maximum  and  the  background.  The  pulse  duration  can  also  be  determined.  However, 
due  to  the  extremely  low  SHG  conversion  efficiency  of  mid-IR  in  nonlinear  crystals, 
the  conventional  setup  is  not  feasible.  To  overcome  this  problem,  we  use  a  two-photon 
quantum  well  infrared  photodetector  (QWIP)  which  converts  the  second-harmonic 
signal  electrically  [36,37]  instead  of  using  a  nonlinear  crystal  plus  a  linear  detector.  The 
energy  diagram  of  one  period  of  the  multi  quantum  well  detector  under  bias  is  shown  in 
the  inset  of  Fig.  11.  The  first  three  electronic  states  are  nearly  equidistant  in  energy. 
When  electrons  in  the  doped  quantum  wells  absorb  two  photons  simultaneously  and  the 
detector  is  biased  (1-3  V),  a  photocurrent  is  generated  and  the  signal  can  be  detected  by 
use  of  a  pre-amplifier  and  a  lock-in  amplifier.  This  experimental  setup  is  diagramed  in 
Fig.  11.  The  second-order  autocorrelation  trace  of  the  MBE-grown  device  mentioned 
above  is  shown  in  Fig.  12.  Interference  fringes  are  observed  when  the  delay  time  from 
one  arm  of  the  autocorrelator  is  equal  to  the  multiples  of  the  cavity  roundtrip  time.  The 
ratio  between  the  maximum  of  the  interference  fringes  and  the  background  is  smaller 


than  8,  and  the  autocorrelation  trace  has  some  features  between  the  cavity  roundtrip 
times,  indicating  that  the  multimode  regime  observed  in  this  device  is  not  stable  mode 
locking  with  a  single  pulse  per  roundtrip. 


12.  A  second-order  autocorrelation  trace  of  a  8  Jim  wavelength  ridge  QC  laser  (wafer  #  2743) 
under  the  condition  of  RNGH  instability.  (Inset)  Microwave  spectrum  of  photocurrent  generated 
by  a  similar  laser  (wafer  #  2721)  under  the  condition  of  RNGH  instability  (measured  with  a  68  kHz 

resolution  bandwidth). 

In  addition  to  second-order  autocorrelation,  the  microwave  spectrum  of  the  laser  output 
was  also  measured  with  an  ultrafast  QWIP  [38]  whose  bandwidth  is  52  GHz.  The  laser 
output  is  sent  directly  to  the  ultrafast  QWIP,  and  the  resulting  photocurrentis  displayed 
in  a  spectrum  analyzer.  Fig.  12  shows  the  beat  note  signal  of  the  MOVPE-grown 
sample  at  pump  current  800  mA  at  77K.  A  steady  peak  with  FWHM  of  13  MHz  at 
22.01  GHz,  which  corresponds  to  the  cavity  round  trip  frequency  of  the  2  mm  long 
laser  (background  refractive  index  n=3),  is  observed  on  the  spectrum  analyzer.  It 
indicates  a  modulation  of  the  laser  output  at  the  cavity  roundtrip  frequency,  and  thus  at 
least  partial  phase  locking  between  the  longitudinal  modes:  The  phase  relationships 
between  the  modes  are  stable  for  about  103  roundtrips. 
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13.  Optical  spectra  vs.  pumping  ratio  ([Trial  mode])  above  threshold  obtained  in  CW  at  77K 
with  a  15  Jim  wide  ridge  laser  emitting  at  8  Jim  (wafer#  2721). 


As  the  ridge  width  is  increased,  the  RNGH  instability  and  the  Rabi  splitting  in  the 
spectrum  are  suppressed.  Fig.  13  shows  the  spectra  of  a  device  processed  from  the  same 
wafer  as  the  one  in  Fig.  9,  with  the  difference  of  its  ridge  width  being  increased  to  15 
pm.  The  spectra  do  not  broaden  much  even  at  very  high  pumping  current,  and  do  not 
show  the  Rabi  splitting.  The  effect  of  Kerr-lensing  decreases  significantly  when  the 
active  region  width  is  increased.  Thus,  this  is  a  strong  evidence  that  the  saturable 
absorption  needed  for  lowering  the  RNGH  threshold  is  provided  by  the  Kerr-lensing 
mechanism. 

5.  Experimental  results:  Buried  heterostructure  lasers 

The  second  type  of  lasers  we  have  tested  are  the  so-called  buried  heterostructure  lasers, 
in  which  an  insulating  Fe-doped  InP  layer  is  regrown  after  etching  of  the  ridges.  A  thick 
layer  of  electrically  plated  gold  is  deposited  on  top  of  the  ridges  as  top  contact  layer 
after  the  InP  regrowth  .  Figure  14  shows  the  cross  section  of  such  a  laser. 


14.  SEM  image  of  the  cross  section  of  a  buried  heterostructure  QC  laser. 

The  active  region  of  the  buried  heterostructure  lasers  tested  is  based  on  the 
four-quantum-well  design,  which  relies  on  a  double  phonon  resonance  to  achieve 
population  inversion  [39].  Fig.  15  a  shows  the  V-I  and  L-I  characteristics  of  a  laser  with 
active  region  width  of  3  pm  and  wavelength  8.38  pm  operated  in  CW  at  room 
temperature,  and  b  shows  its  optical  spectra.  The  spectra  show  Rabi  splitting  similar  to 
the  ridge  laser  case,  indicating  the  RNGH  instability  in  this  narrow  buried 
heterostructure  laser.  The  spectral  splitting  and  twice  the  Rabi  frequency  are  plotted 
against  the  square  root  of  the  collected  output  power  in  Fig.  15  b.  A  good  agreement  is 
found  between  the  experimental  splitting  and  twice  the  calculated  Rabi  frequency. 


so 


FIG.  15:  (a)  V-I  &  L-I  curves  and  (b)  optical  spectra  vs.  pumping  ratio  {j/jth)  above  threshold  obtained  in  CW  at  300K  with 
a  3  pm  wide  ridge  laser  emitting  at  8.38  pm  (wafer  #  3251).  (c)Spectral  splitting  and  twice  the  Rabi  frequency  fiRabt/(27r)  vs. 
square  root  of  output  power  collected  from  a  single  laser  facet.  The  dashed  line  is  a  least-square  linear  fit  of  the  data. 

The  second-order  autocorrelation  trace  of  the  device  when  it  is  pumped  five  times 
above  threshold  at  80K  is  shown  in  Fig.  16.  The  ratio  between  the  maximum  of  the 
interference  fringes  and  the  background  is  close  to  8  to  3,  similar  to  the  ridge  laser  case. 
In  addition,  there  are  smaller  interference  fringes  within  one  cavity  roundtrip  period. 
Both  features  indicate  no  stable  pulsation  from  the  laser.  The  microwave  spectrum  of 
the  laser  output  shows  a  steady  peak  with  FWHM  of  700  kHz  at  15.018  GHz.  This  is 
shown  in  the  inset  of  Fig.  16. 


Delay  [psj 


16.  A  second-order  autocorrelation  trace  of  a  8  Jim  wavelength  buried  heterostructure  QC  laser 
(wafer  #  3251)  under  the  condition  of  RNGH  instability.  (Inset)  Microwave  spectrum  of  the 
photocurrent  generated  by  the  same  laser  (measured  with  a  68  kHz  resolution  bandwidth). 

As  in  the  ridge  laser  case,  the  lowering  of  the  threshold  of  RNGH  instability  can  also  be 
attributed  to  the  Kerr-lens  type  saturable  absorption.  In  order  to  better  support  this  idea 


experimentally,  the  spectra  from  another  device  processed  from  the  same  wafer  as  in 
Fig.  15  but  with  a  wider  active  region  (7.5  pm)  are  measured.  Two-dimensional 
waveguide  simulations  indicate  a  much  weaker  Kerr-lensing  effect  in  these  QCLs,  due 
to  the  much  larger  ratio  of  active  region  width  to  wavelength.  The  measured  optical 
spectra  obtained  at  300K  in  CW  mode  are  shown  in  Fig.  17  a.  The  envelopes  of  the 
spectra  consist  of  multiple  peaks  whose  separation  is  independent  of  the  pumping 
current.  The  spectral  signatures  qualitatively  agree  with  the  numerical  simulations  in  a 
Fabry-Perot  cavity  without  a  saturable  absorber  (y  =0)  (Fig.  6). 
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FIG.  17:  Optical  spectra  vs.  pumping  ratio  (j/jth)  above  threshold  obtained  in  CW  at  300K  with  (a)  a  7.5  pm  wide  and 
a  10  pm  buried  heterostructure  lasers  emitting  at  8  pm  (wafer  #  3251) 


Further  increasing  the  active  region  width  tends  to  suppress  both  the  instabilities  caused 
by  RNGF1  and  spatial  hole  burning.  Fig.  17  b  shows  the  spectra  of  a  device  also 
processed  from  the  same  wafer  but  with  an  even  wider  active  region  (10  pm)  at  77K. 
The  spectra  do  not  broaden  much  even  at  very  high  pumping  currents,  as  in  the  ridge 
laser  case.  Since  both  RNGF1  instability  and  SF1B  stem  from  nonlinear  effects,  they  are 
suppressed  when  the  intensity  of  the  field  in  the  cavity  is  lower.  Moreover,  when  the 
active  region  width  increases,  higher  transverse  modes  are  excited.  Different  transverse 
modes  have  different  propagation  constants  [3,  and  thus  form  different  gain  gratings 
which  tend  to  wash  out  the  effect  of  spatial  hole  burning. 


6.  Temperature  effects 


FIG.  18:  Optical  spectra  vs.  pumping  ratio  (j/jth)  above  threshold  obtained  in  cw  of  a  buried  heterostructure  laser  emitting 
at  8.38  pm  at  different  temperatures:  (a)80K,  (b)150K,  (c)200K.  and  (d)300K 


Temperature  also  plays  an  important  role  in  the  shape  of  the  spectrum.  To  illustrate  this 
point,  we  now  present  the  spectra  at  different  temperatures  for  the  a=8.38  pm  buried 
heterostructure  laser  with  3  pm  active  region  width  (the  same  as  in  Fig.  15). 

From  the  spectra  in  Fig.  18,  it  is  clear  that  at  lower  temperatures  the  spectra  are 
dominated  by  spatial  hole  burning,  showing  multiple  peaks  independent  of  pumping 
and  no  significant  Rabi  splitting.  As  the  temperature  increases,  the  Rabi  splitting 
becomes  more  evident  and  finally  the  RNGFI  instability  takes  over. 

Intuitively,  the  effect  of  the  temperature  on  the  nature  of  the  multimode  regime  in  QCLs 
seems  to  come  from  carrier  diffusion.  The  lifetime  of  the  gain  grating  Tg  is  defined  as 
Tg  1  =  T~l  +  4k2 D .  The  diffusion  coefficient  D  of  the  gain  grating  is  proportional  to  the 

temperature,  D  =  fjkBT / q  ,  where  p  here  is  the  carrier  mobility,  ke  is  the  Boltzmann’s 
constant,  T  is  the  temperature  and  q  is  the  carrier  charge.  Therefore  at  higher 
temperature  carrier  diffusion  would  reduce  spatial  hole  burning  and  thus  reveal  the 
RNGFI  instability.  However,  with  a  mobility  of  16000  cm2/sec/V  and  k=2.25*104  cm"1 
(which  corresonds  to  a  vacuum  wavelength  of  8.38pm),  4k2D  ~  0.2  THz  at  77K  and 
4k2D  ~  0.8  THz  at  300K,  both  significantly  smaller  than  T f1  ~  0.6  THz  at  77K  and  T f1 
~  2  THz  at  300K.  Thus,  carrier  diffusion  is  unlikely  to  be  the  reason  for  the  temperature 
dependence.  Although  the  temperature  effect  is  not  entirely  understood  to  us  at  this 
point,  one  possiblity  is  temperature-dependent  saturable  absorption.  The  QCL  injector 
consists  of  many  energy  levels  which  can  be  thermally  populated.  It  is  not  surprising 
that  if  any  two  higher  levels  in  the  injector  are  closely  resonant  with  the  laser  transition, 
and  that  they  will  form  a  two-level  saturable  absorber.  Therefore  in  this  case  saturable 
absorption  is  stronger  at  higher  temperature  and  makes  the  RNGH  more  easily 


observable. 


7.  Conclusion 

This  paper  provides  a  thorough  account  of  different  multimode  regimes  in  QCLs.  It  was 
found  that  two  key  mechanisms  which  govern  the  multimode  regimes  in  QCLs  are  a 
coherent  instability  similar  to  the  RNGH  instability,  and  spatial  hole  burning.  Both 
mechanisms  are  enhanced  due  to  the  large  dipole  moment  p  of  the  laser  transition, 
which  results  in  the  unusually  fast  gain  recovery  in  QCLs:  The  RNGH  instability  is 
enhanced  because  the  Rabi  frequency  scales  as  p,  and  therefore  the  Rabi  splitting  can 
be  resolved  by  the  comb  of  modes  supported  by  the  cavity.  SHB  is  enhanced  because 
carrier  diffusion  is  slower  than  the  gain  recovery,  and  thus  leaves  the  gain  grating 
intact. 

Due  to  the  fast  gain  recovery,  conventional  mode  locking,  with  one  pulse  per  roundtrip, 
is  suppressed.  In  order  to  achieve  conventional  mode  locking  in  QCLs,  one  needs  to 
design  a  QCL  with  a  slower  gain  recovery,  such  that  Tj  becomes  longer  than  or 
comparable  to  the  cavity  roundtrip.  Efforts  in  this  direction  are  currently  underway. 
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A  theoretical  and  experimental  study  of  multimode  operation  regimes  in  quantum  cascade  lasers  (QCLs)  is 
presented.  It  is  shown  that  the  fast  gain  recovery  of  QCLs  promotes  two  multimode  regimes:  One  is  spatial 
hole  burning  (SHB)  and  the  other  one  is  related  to  the  Risken-Nummedal-Graham-Haken  instability  predicted 
in  the  1960s.  A  model  that  can  account  for  coherent  phenomena,  a  saturable  absorber,  and  SHB  is  developed 
and  studied  in  detail  both  analytically  and  numerically.  A  wide  variety  of  experimental  data  on  multimode 
regimes  is  presented.  Lasers  with  a  narrow  active  region  and/or  with  metal  coating  on  the  sides  tend  to  develop 
a  splitting  in  the  spectrum,  approximately  equal  to  twice  the  Rabi  frequency.  It  is  proposed  that  this  behavior 
stems  from  the  presence  of  a  saturable  absorber,  which  can  result  from  a  Kerr  lensing  effect  in  the  cavity. 

Lasers  with  a  wide  active  region,  which  have  a  weaker  saturable  absorber,  do  not  exhibit  a  Rabi  splitting  and 
their  multimode  regime  is  governed  by  SHB.  This  experimental  phenomenology  is  well-explained  by  our 
theoretical  model.  The  temperature  dependence  of  the  multimode  regime  is  also  presented. 
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I.  INTRODUCTION 

Many  physical  mechanisms  can  drive  a  laser  from  a 
single-mode  to  a  multimode  regime.  Common  examples  are 
spatial  and  spectral  hole  burning  (SHB),  saturable  absorp¬ 
tion,  and  self-phase  modulation  [1-3].  Understanding  these 
mechanisms  is  of  key  importance  to  laser  science  and  tech¬ 
nology,  whether  one  is  interested  in  single-mode  behavior  or 
in  a  particular  multimode  operation  such  as  mode  locking. 

The  multimode  regimes  listed  above  are  well-understood 
and  documented  both  theoretically  and  experimentally  in 
many  types  of  lasers.  However,  the  understanding  of  multi- 
mode  regimes  in  quantum  cascade  lasers  (QCLs)  is  still  in  its 
infancy,  as  these  lasers  were  only  demonstrated  in  1994  [4], 
and  studies  of  their  multimode  regimes  commenced  even 
more  recently  [5-8], 

As  it  was  shown  recently  [9]  multimode  dynamics  in 
QCLs  is  different  from  that  of  more  common  lasers.  This  is 
mainly  due  to  the  unusually  fast  gain  recovery  of  QCLs, 
which  occurs  on  a  picosecond  scale.  While  a  saturable  ab¬ 
sorber  triggers  mode  locking  in  lasers  with  slow  gain  recov¬ 
ery  (relative  to  the  round-trip  time),  in  lasers  with  fast  gain 
recovery  a  saturable  absorber  triggers  a  mechanism  [9]  simi¬ 
lar  to  the  Risken-Nummedal-Graham-Haken  (RNGH)  insta¬ 
bility  [10,11],  While  in  standard  semiconductor  lasers  carrier 
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diffusion  eliminates  spatial  hole  burning,  in  QCLs  the  gain 
recovery  process  is  faster  than  carrier  diffusion,  and  spatial 
hole  burning  is  dominant,  favoring  multimode  operation. 

This  paper  presents  a  detailed  study  of  multimode  regimes 
in  QCLs.  The  results  of  Ref.  [10,11]  are  substantiated  and 
extended.  The  first  part  of  the  paper  is  theoretical  and  the 
second  one  is  experimental.  In  the  theoretical  section,  the 
laser  Maxwell-Bloch  equations  are  introduced  and  analyzed 
for  both  a  ring  laser  cavity  and  a  Fabry-Perot  cavity.  This 
way  one  can  study  the  interplay  of  coherent  phenomena  and 
spatial  hole  burning.  A  saturable  absorber  is  added  to  the 
model  as  well.  The  model  is  studied  analytically  and  numeri¬ 
cally.  The  stability  region  of  a  continuous  wave  (cw)  solution 
is  found.  It  is  shown  that  in  a  ring  cavity,  the  presence  of  a 
saturable  absorber  lowers  the  threshold  of  the  RNGH  insta¬ 
bility  from  about  nine  times  above  laser  threshold  [10,11]  to 
arbitrarily  low  above  laser  threshold,  depending  on  the 
strength  of  the  absorber.  However,  the  nature  of  the  instabil¬ 
ity  remains  the  same:  the  population  inversion  begins  to  os¬ 
cillate  at  the  Rabi  frequency,  modulating  the  gain  in  the  laser. 
The  result  is  sidebands  around  the  original  cw  mode,  sepa¬ 
rated  from  it  by  roughly  the  Rabi  frequency.  A  Rabi  splitting 
in  the  spectrum  is  the  primary  signature  of  the  RNGH  insta¬ 
bility. 

It  is  then  shown  that  the  ring-cavity  Maxwell-Bloch 
model  with  a  saturable  absorber  can  explain  the  Rabi  split¬ 
ting  and  the  lowering  of  the  threshold,  but  cannot  explain  a 
key  feature  in  the  experimental  spectra.  The  latter  exhibit 
two  relatively  equal  groups  of  modes  with  a  gap  in  between. 
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whereas  the  ring-cavity  Maxwell-Bloch  model  predicts  a 
large  central  mode  in  the  spectrum  with  two  sidebands.  In 
view  of  this  discrepancy,  the  Maxwell-Bloch  equations  were 
extended  to  include  coupling  between  counterpropagating 
modes  in  a  Fabry-Perot  cavity,  which  is  of  course  a  more 
adequate  model  for  QCLs.  This  model  takes  into  account  the 
development  of  SHB.  The  inclusion  of  SHB  generates  theo¬ 
retical  spectra  with  a  Rabi  splitting  and  without  a  central 
peak  in  the  spectrum,  in  agreement  with  experiments. 

After  presenting  the  theoretical  results,  a  thorough  study 
of  the  experimental  phenomenology  is  presented.  It  is  shown 
that  QCLs  with  narrower  active  regions  tend  to  exhibit  a 
more  pronounced  Rabi  splitting  than  lasers  with  wider  active 
regions.  Lasers  with  a  wider  active  region  tend  to  exhibit 
multimode  spectra  that  are  governed  by  SHB.  The  explana¬ 
tion  we  propose  to  this  behavior  is  that  narrow  QCLs  have  a 
stronger  saturable  absorption  effect  than  wider  ones  because 
the  optical  intensity  (and  thus  any  nonlinear  effect)  is  en¬ 
hanced,  and  as  a  result  of  an  enhanced  Kerr-lensing  effect 
[3,5].  In  this  case  a  nonlinear  index  enhancement  in  the 
waveguide  core  gives  rise  to  an  increased  overlap  of  the 
transverse  laser  mode  with  an  active  region  and  a  reduced 
overlap  with  lossy  sidewalls,  leading  to  an  additional  en¬ 
hancement  of  the  saturable  absorption.  The  idea  of  saturable 
absorption  by  Kerr  lensing  is  further  supported  by  the  fact 
that  QCLs  with  metal  coating  on  the  sides  of  the  ridge  have 
a  stronger  RNGH  behavior  than  lasers  without  metal  coating: 
The  metal  coating  enhances  saturable  losses  that  originated 
from  Kerr  lensing. 

In  the  last  part  of  the  paper,  the  temperature  dependence 
of  the  multimode  behavior  in  QCLs  is  studied  experimen¬ 
tally.  It  is  found  that  at  higher  temperature  the  multimode 
behavior  tends  to  be  governed  by  the  RNGH  instability, 
whereas  at  lower  temperatures  it  is  governed  by  SHB.  This 
behavior  may  be  explained  by  the  fact  that  at  higher  tem¬ 
peratures  hot  carriers  populate  more  states  in  the  injector 
superlattice,  creating  additional  quasiresonant  absorption 
transitions  between  ground  and  excited  minibands.  This 
leads  to  an  additional  saturable  absorption  of  laser  radiation. 

This  paper  is  organized  as  follows:  Section  II  gives  a  brief 
survey  of  prior  work  on  multimode  regimes  in  QCLs,  Secs. 
III-V  study  theoretically  the  Maxwell-Bloch  equations  in  a 
ring  cavity  and  a  Fabry-Perot  cavity.  Secs.  VI-VIII  summa¬ 
rize  the  experimental  study,  and  Sec.  IX  is  a  brief  summary. 


II.  PRIOR  WORK  ON  MULTIMODE  REGIMES 
IN  QCLS 

Multimode  regimes  in  QCLs  were  observed  in  a  series  of 
recent  works  [5-8].  In  Ref.  [5],  for  example,  it  was  observed 
that  at  a  certain  pumping  current  above  lasing  threshold, 
QCLs  cease  to  operate  in  cw  and  develop  a  multimode  re¬ 
gime.  This  multimode  regime  was  characterized  by  a  broad¬ 
band  optical  spectrum  and  a  narrow  (less  than  100  kHz)  ra¬ 
dio  frequency  (rf)  beat  note  in  the  power  spectrum  at  the 
cavity  round-trip  frequency. 

The  narrow  beat  note,  whose  width  is  1  / 105  of  the  central 
frequency,  shows  that  the  wave  form  of  the  electric  field 
circulating  in  the  laser  cavity  was  stable  over  approximately 


105  round-trips.  In  other  words,  the  phase  relationships  be¬ 
tween  the  longitudinal  modes  were  stable  for  about  105 
round-trips.  The  modes  were  therefore  locked.  However,  in 
order  to  characterize  the  wave  form  which  is  circulating  in 
the  laser  and  to  see  if  it  is  indeed  an  isolated  pulse,  as  in 
traditional  mode  locking,  one  has  to  apply  pulse  character¬ 
ization  techniques,  such  as  second  order  autocorrelation. 

At  the  time  when  the  experiment  in  Ref.  [5]  was  per¬ 
formed,  no  second-order  autocorrelation  apparatus  was  avail¬ 
able.  However,  second  harmonic  generation  from  QCLs  [6] 
provided  some  information  for  pulse  characterization.  This 
measurement  showed  an  increase  by  more  than  a  factor  of  5 
in  the  second  harmonic  signal  as  the  multimode  behavior  set 
in.  This  increase  indicates  that  the  duty  cycle  of  the  pulses 
was  roughly  5.  Since  the  number  of  modes  in  the  spectrum  is 
more  than  5,  we  conclude  that  not  all  the  modes  were  locked 
into  a  pulse.  Without  better  pulse  characterization  data,  one 
can  infer  that  the  laser  could  have  had  more  than  one  pulse 
per  round-trip. 

It  should  also  be  noted  that  traditional  mode  locking,  with 
a  single  pulse  per  round-trip,  requires  that  the  gain  recovery 
time  be  longer  than  the  cavity  round-trip  [3].  In  QCLs  this 
condition  is  violated,  and  therefore  according  to  the  existing 
theory  one  cannot  expect  mode  locking  with  a  single  pulse 
per  round-trip  in  QCLs.  One  can  expect  multiple  pulses  per 
round-trip. 

In  view  of  the  above,  one  can  see  that  the  nature  of  vari¬ 
ous  multimode  regimes  in  QCLs  requires  further  elucidation. 
This  is  the  main  goal  of  the  present  work. 


III.  MAXWELL-BLOCH  EQUATIONS  IN  A  FABRY-PEROT 
CAVITY 

In  this  section  we  derive  the  Maxwell-Bloch  equations  in 
a  Fabry-Perot  cavity.  We  model  the  gain  medium  of  QCLs  as 
a  two  level  system,  described  by  the  Bloch  equations 

■  .  Pab  /. .. 

Pab  =  'Mpab+l  —  !\-—,  (1) 

ft  T2 


2i^~(p"h  -  Pab)  ■ 


A„  -  A 


+  D~ 


oPA 

*2’ 
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where  pab  is  the  off-diagonal  element  of  the  density  matrix, 
A  =  Pi,h  ~  Pan  is  the  population  inversion,  o>  and  p  are  the 
resonant  frequency  and  the  dipole  matrix  element  of  the  las¬ 
ing  transition,  7)  and  T2  are  the  longitudinal  and  transverse 
relaxation  times,  and  Ap  is  the  steady-state  inversion  at  E 
=  0,  which  characterizes  the  pumping  rate.  The  last  term  in 
Eq.  (2)  is  added  phenomenologically  and  accounts  for  spatial 
diffusion  of  the  inversion  due  to  carrier  diffusion.  D  is  the 
diffusion  coefficient.  E  is  the  electric  field,  which  is  assumed 
to  satisfy  the  wave  equation 

dLzE-n-2dL,E=^dlt{pab  +  p*ab).  (3) 

c  e0c 


N  is  the  number  of  two-level  systems  per  unit  volume,  which 
equals  the  average  doping  density  in  the  active  regions.  F  is 
the  overlap  factor  between  the  optical  mode  and  the  active 
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region,  and  n  is  the  background  refractive  index. 

Equations  (1)  and  (2)  describe  an  open  two-level  system 
[12]:  The  total  number  of  electrons  in  the  system  is  not  con¬ 
served,  but  rather  can  flow  in  and  out,  and  vary  depending  on 
the  bias  conditions.  In  general  the  two  levels,  the  upper  one 
b  and  the  lower  one  a ,  can  have  different  relaxation  times.  In 
this  situation  Eq.  (2)  should  be  replaced  by  two  equations, 
each  with  its  own  relaxation  times.  In  QCLs  indeed  the  upper 
and  lower  levels  have  different  relaxation  times.  However, 
we  neglect  this  difference  for  simplicity’s  sake  and  assign  to 
them  the  same  value  7],  The  generalization  of  the  model  (1) 
and  (2)  to  a  model  with  two  different  relaxation  times  is 
straightforward,  as  is  also  the  generalization  to  a  model  with 
more  than  two  levels. 

We  now  make  the  following  set  of  ansatzs: 

E(z,t)  =  ^[El(z,t)e-{ia'-kz)  +  £+(z,f)e/(“'-te)] 

+  {-\E*Sz,i)e-{iM+kz)  +  E_(z,t)ei{M+kz)],  (4) 

pab(z,  t )  =  77 +(z,  t)ei{w,-kz)  +  77_(z,  i)ei(a,t+kz) ,  (5) 

A  (z,t)  =  A0(z,f)  +  A  2(z,t)e2ikz  +  A  *2(z,t)e~2ikz,  (6) 

where  k=amlc.  E±,  t]±,  A0,  and  A2  are  assumed  to  vary 
slowly  in  time  and  space,  on  the  scale  defined  by  a>  and  k, 
respectively.  The  quantities  with  a  +  (-)  subscript  represent 
waves  traveling  to  the  positive  (negative)  z  direction.  Equa¬ 
tion  (6)  allows  taking  SHB  into  account,  A2  being  the  enve¬ 
lope  of  the  inversion  grating.  Note  that  Eq.  (6)  can  be  ex¬ 
tended  by  adding  terms  proportional  to  eAlkz,  e6,kz,  etc. 
Neglecting  these  terms  means  that  higher  spatial  frequencies 
on  the  inversion  are  neglected.  Due  to  the  nonlinearity  of  the 
gain  medium  such  frequencies  can  appear  if  the  gain  is 
heavily  saturated,  but  they  are  neglected  in  our  analysis  for 
simplicity’s  sake  and  since  the  pumping  in  our  system  is 
never  much  higher  than  the  laser  threshold. 

Substituting  Eqs.  (4)-(6)  into  Eqs.  (l)-(3)  and  making  the 
slowly  varying  envelope  approximation,  we  obtain  the  fol¬ 
lowing  set  of  equations: 

n  kNfiT  1 

-d,E±  =  +  dzE±  -  V±  -  -£{E+,E_)E±,  (7) 

c  2  e0n~  2 

dlrf±='^(AQE±  +  AfET)-^,  (8) 

2  n  T2 

^Ao  =  +  l~~{E*+r]+  +  E*_r/_  -  c.c.),  (9) 

T\  n 

dAt  =  ±  ijiKv-  ~  V+E-)  -  ^  -  4k2  DA2 ,  (10) 

n  T{ 

where  we  have  introduced  the  notation  A2  =  A2,  A2  =  A2,  in 
order  that  Eqs.  (7)-(10)  can  be  written  more  compactly.  The 
last  term  in  Eq.  (7)  has  been  added,  and  represents  loss.  The 
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loss  l  is  allowed  to  depend  on  the  field  to  represent  phenom¬ 
ena  such  as  optical  saturation. 

In  QCLs  the  laser  cavity  is  formed  by  the  two  cleaved 
facets,  one  located  at  z=0  and  the  other  one  at  z=L.  At  each 
facet  the  Fresnel  reflection  law  dictates  the  following  rela¬ 
tions: 


E_{L,t)-  E+{L,t), 

77+1 

(11) 

n  —  1 

EM-  EM- 

77+1 

(12) 

In  what  follows  we  study  analytically  and  numerically  the 
model  introduced  in  this  section.  We  begin  with  a  simplified 
case,  namely  the  standard  Maxwell-Bloch  equation  in  a  ring 
cavity.  In  addition  to  briefly  reviewing  known  results  about 
the  RNGH  instability,  we  study  the  effect  of  a  saturable  ab¬ 
sorber  on  the  latter. 


IV.  RING  CAVITY 

In  this  section  we  consider  a  ring  cavity,  where  SHB  does 
not  exist  because  standing  waves  cannot  form.  The  aim  is  to 
understand  the  interplay  of  coherent  effects  and  a  saturable 
absorber  alone,  while  avoiding  complications  due  to  SHB. 
We  shall  see  that  without  SHB,  the  qualitative  agreement 
between  theory  and  the  experiments  on  QCLs  is  not  com¬ 
plete.  After  introducing  SHB  in  the  next  section,  the  agree¬ 
ment  is  much  more  satisfactory. 


A.  RNGH  instability  with  a  saturable  absorber 

Dropping  all  the  quantities  with  a  subscript  from  Eqs. 
(7)-(10),  one  arrives  at  the  standard  Maxwell-Bloch  equa¬ 
tions,  with  a  saturable  absorber  added. 

~d,E  =  -  dzE-^--\{lQ-1\E\2)E,  (13) 

c  hi  Q  Ath  2 


d,rj-- 


ILL  71 

— AE-— , 
2  h  T2 


(14) 


9,  A  = 


Tx 


(E*7]~  C.C.). 


(15) 


The  saturable  absorber  is  approximated  to  the  lowest  order  in 
E,  and  is  characterized  by  y,  which  is  often  referred  to  as  the 
self-amplitude  modulation  coefficient  [3].  l0  is  the  linear 
loss.  Ath  is  the  lasing  threshold  value  of  Ap  for  y=0,  given  by 
(see  Appendix  A) 


■i  _  kN]i2TT2 
1  2M0e0n2 


(16) 


Linear  stability  analysis  of  Eqs.  (13)— (15)  is  carried  out  in 
detail  in  Appendix  A.  The  gain  of  a  perturbation  at  the  fre¬ 
quency  11  (relative  to  the  resonance  frequency  «)  is  approxi¬ 
mately  given  by 
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FIG.  1.  (Color  online)  g(fl)  for  p= 2  and  the  parameters  in  Table 
I,  apart  from  y.  The  latter  is  y=  0  (solid),  y=10-9  cm/V2  (dashed), 
and  y=2X  10-9  cm/V2  (dotted). 


=  -  —Re 
2  n 


(fir1  +  /)nr2-2(^-i) 

\(lT1+i){D,T2  +  i)-(p-  1) 


yh\p  - 1)  (nr,  +  0(3flr2  +  2/)  -  4(P  - 1) 
p}TxT2  (ilT |  +  i)(ilT2  +  i)  - p  +  1 


(17) 


The  approximations  made  in  the  derivation  of  Eq.  (17)  are 
discussed  in  detail  in  Appendix  A  and  mainly  include  assum¬ 
ing  that  the  photon  lifetime  in  the  empty  cavity  is  much 
longer  than  Tx  and  72.  This  approximation  is  excellent  for 
QCLs.  p  is  the  pumping  above  lasing  threshold  (for  y=0): 


Figure  1  shows  #(0)  for  different  values  of  y  at  p= 2.  The 
parameters  used  in  all  examples  in  this  paper,  unless  speci¬ 
fied  otherwise,  are  given  in  Table  I.  The  effect  of  y  is  to 
increase  g(fl)  more  or  less  uniformly  across  the  frequency 
domain.  In  particular,  it  can  bring  g(H)  above  zero,  thereby 
triggering  an  instability,  even  when  the  laser  is  not  pumped 
as  high  above  threshold.  The  reason  why  a  saturable  absorber 
lowers  the  RNGH  threshold  is  that  a  saturable  absorber  itself 
always  favors  a  multimode  regime  to  a  single  mode  one.  It 
introduces  a  frequency-independent  parametric  gain.  The  lat¬ 
ter  is  added  to  the  RNGH  parametric  gain  from  Fig.  1,  bring¬ 
ing  it  above  threshold.  Note,  however,  that  the  instability  still 


TABLE  I.  The  parameters  used  in  all  calculations  and  simula¬ 
tions  in  this  paper,  unless  indicated  otherwise. 


Quantity 

Symbol 

Value 

Gain  recovery  time 

Ty 

0.5  ps 

Dephasing  time 

t2 

0.067  ps 

Linear  cavity  loss 

^0 

5  cm-1 

Transition  dipole  element 

A 

2.54  nmXe 

Background  refractive  index 

n 

3.3 

Cavity  length 

L 

3  mm 

Saturable  absorber  coefficient 

r 

10“8  cm/V2 

FIG.  2.  (Color  online)  The  pumping  ratio  p  at  which  the  RNGH 
instability  sets  in  as  a  function  of  the  saturable  absorber  coefficient 
y.  Solid  line:  exact  result  from  linear  stability  analysis.  Dashed  line: 
approximate  result  [Eq.  (21)].  The  parameters,  apart  from  y,  are 
given  in  Table  I. 

starts  by  developing  Rabi  sidebands  around  the  cw  lasing 
frequency.  In  this  sense,  it  can  be  interpreted  as  a  modified 
version  of  the  original  RNGH  instability,  rather  than  as  a 
modulation  instability  caused  by  the  saturable  absorber 
alone. 

The  dependence  of  the  threshold  for  instability  on  y  is 
shown  in  Fig.  2.  For  y=0  we  recover  the  standard  RNGH 
instability,  which  occurs  when  pumping  of  slightly  above  9 
times  the  lasing  threshold.  For  y— >°°,  the  instability  thresh¬ 
old  approaches  the  lasing  threshold. 

The  results  in  Figs.  1  and  2  were  obtained  from  numerical 
solutions  of  the  algebraic  equations  involved  in  the  stability 
analysis  (see  Appendix  A).  Since  even  the  approximate  ex¬ 
pression  of  g(fl)  [Eq.  (17)]  is  not  very  simple,  it  is  useful  to 
derive  some  approximate  simple  expressions  for  the  proper¬ 
ties  of  the  instability.  This  is  done  in  detail  in  Appendix  A, 
and  here  we  only  give  the  results. 

g(fl)  (e.g.,  in  Fig.  1)  has  a  local  minimum  at  fl=0,  and 
peaks  at  approximately 

l^maxl  =  ftRab .  V  ^7  ■  (19) 

V  p  -  1 

The  Rabi  frequency  0R.lbl  =  piElh,  where  pc  is  the  electron 
charge  times  the  matrix  element  of  the  laser  transition.  The 
position  of  the  peak  depends  weakly  on  y  within  the  param¬ 
eter  range  of  interest  to  our  system.  Note  that  for  y=  0  one 
has 


l^Rabi  - 


T\T2 


(20) 


The  dependence  of  |Hmax|  on  p  is  shown  in  Fig.  3. 
The  instability  threshold  is  approximately  given  by 


Pth  -  1  +  8 


h2y 


2  h2y 

+  12^W 

p,  1  2to 


+  1 


(21) 


Figure  2  shows  that  Eq.  (21)  fairly  well  approximates  the 
exact  threshold  condition. 
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FIG.  3.  (Color  online)  Hmax  as  a  function  of  the  pumping.  Thick 
line:  exact,  and  dotted  line:  Eq.  (19).  The  parameters  are  from  Table 
I  with  the  only  difference  that  y=0. 

B.  Numerical  results 

The  linear  stability  analysis  given  above  can  only  find  the 
threshold  condition,  and  does  not  describe  the  behavior  of 
the  laser  above  the  point  of  instability.  Here  we  therefore 
present  the  results  of  a  numerical  solution  of  Eqs.  (13)— (15). 

According  to  Ref.  [13],  the  energy  in  the  Rabi  sidebands 
can  change  either  discontinuously  or  continuously  at  the 
RNGH  instability  threshold.  In  the  language  of  phase  transi¬ 
tions,  this  would  be  a  first-  or  second-order  phase  transition, 
respectively.  The  criterion  that  determines  the  order  of  the 
transition  involves  the  laser  cavity  length,  but  has  not  been 
given  explicitly  in  the  literature  known  to  us. 

According  to  the  numerical  simulations,  for  the  param¬ 
eters  corresponding  to  our  QCLs,  including  the  saturable  ab¬ 
sorber,  the  transition  is  continuous  (second  order).  Below  the 
instability  threshold  the  lasing  is  cw,  and  as  the  threshold  is 
crossed,  the  Rabi  sidebands  around  the  central  cw  mode 
grow  continuously.  This  behavior  is  demonstrated  in  Fig.  4. 
Note  that  standard  mode  locking,  with  a  saturable  absorber 
and  slow  gain  recovery,  is  a  discontinuous  (first-order)  tran¬ 
sition  [14,15]. 


The  spectra  in  Fig.  4  have  two  groups  of  modes  separated 
by  roughly  twice  the  Rabi  frequency.  In  this  respect  they 
resemble  the  experimental  spectra.  However,  they  have  a 
strong  cw  peak  in  between,  a  feature  which  is  not  shared  by 
the  experimental  spectra.  As  we  show  in  the  next  section,  in 
a  Fabry-Perot  cavity  with  SHB,  the  central  cw  peak  disap¬ 
pears. 

V.  FABRY-PEROT  CAVITY 

Spatial  hole  burning  is  associated  with  A2  in  Eq.  (10). 
Intuitively,  A2  is  the  amplitude  of  the  grating  that  couples  the 
two  propagation  directions  in  the  laser.  The  parameter  that 
controls  the  strength  of  SHB  is  D:  in  the  limit  of  D— »°°,  A2 
approaches  zero.  In  order  to  better  understand  the  interplay 
between  SHB  and  the  RNGH  instability,  we  present  now  the 
results  of  analytical  and  numerical  studies  of  Eqs.  (7)-(10). 
We  start  with  linear  stability  analysis.  The  calculation  is 
shown  in  detail  in  Appendix  A,  and  here  we  only  give  the 
results. 

Before  proceeding,  from  Eq.  (10)  we  define  the  lifetime 
of  the  gain  grating  Tg  as 

T~g  =  T\l  +  4k2D.  (22) 

Tg  is  the  parameter  that  determines  the  strength  of  spatial 
hole  burning.  Tg  can  therefore  range  from  zero  (no  SHB)  to 
7)  (strongest  of  SHB).  The  diffusion  coefficient  D  can  be 
estimated  from  the  Einstein  relation.  With  an  electron  mobil¬ 
ity  of  7000  (cm2/s)/V  at  room  temperature,  one  has  D 
=  180  cm2/s,  k=3.7  X  104  cm  1 ,  which  roughly  corresponds 
to  a  vacuum  wavelength  of  5  fin l,  and  we  obtain  4 k2D 
~1  THz.  With  T\l~2  THz  we  find  that  Tg~ 0.3  ps.  Note 
that  the  mobility  used  here  was  relatively  high,  and  the 
wavelength  was  on  the  short  side  of  the  scale.  Therefore  in 
reality  Tg  is  closer  to  7).  It  therefore  follows  that  due  to  the 
fast  gain  recovery  of  QCLs,  carrier  diffusion  does  not  elimi¬ 
nate  spatial  hole  burning.  This  is  in  contrast  with  diode  lasers 
[16,17]. 

In  single  mode  operation,  the  standing  wave  associated 
with  the  lasing  mode  imprints  a  grating  in  the  medium  which 
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FIG.  4.  (Color  online)  Spectra 
obtained  from  a  numerical  solu¬ 
tion  of  the  Maxwell-Bloch  equa¬ 
tions  [Eqs.  (13)— (15)].  The  param¬ 
eters  are  given  in  Table  I,  with  the 
following  exceptions:  y=  1 

X  10-9  cm/V2,  L= 6  mm  (to  give 
the  same  round-trip  time  as  in  a 
Fabry-Perot  3  mm  long  cavity), 
and  f0=7  cm-1  (to  account  for 
the  mirror  losses  too). 
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FIG.  5.  (Color  online)  g(fl)  of  the  perturbation,  associated  with 
spatial  hole  burning,  p  and  Tg  corresponding  to  every  curve  are 
denoted  in  the  legend. 

has,  to  first  approximation,  a  sinusoidal  profile.  This  grating 
causes  the  lasing  mode  to  experience  heavier  gain  saturation 
than  any  other  mode  would  experience.  When  the  gain  grat¬ 
ing  becomes  strong  enough,  the  single  mode  regime  becomes 
unfavorable,  and  additional  modes  are  excited. 


Optical  Frequency  (THz) 


1900  1925  1950  1975  2000 


Wave  Number  (cm-1) 

FIG.  6.  Results  of  numerical  simulations  of  the  spectra  based  on 
the  Maxwell-Bloch  equations  including  a  saturable  absorber  and 
spatial  hole  burning  for  different  values  of  the  current  density  nor¬ 
malized  to  the  threshold  value. 


A.  Linear  stability  analysis 

Linear  stability  analysis  of  Eqs.  (7)-(10)  gives  two  fami¬ 
lies  of  unstable  modes.  One  is  associated  with  the  RNGH 
instability,  and  the  other  one  with  spatial  hole  burning.  The 
first  family  is  very  similar  to  the  case  of  a  ring  cavity  studied 
earlier,  with  small  differences  that  are  discussed  in  Appendix 
A.  The  second  one  is  derived  in  Appendix  A,  and  the  result  is 
shown  in  Fig.  5. 

Figure  5  shows  the  gain  of  a  perturbation  around  a  cw 
solution.  The  shorter  is  Tg  and  the  smaller  is  the  pumping, 
the  smaller  is  the  gain  of  the  instability.  For  parameters  typi¬ 
cal  to  our  QCLs,  the  SHB  instability  occurs  a  few  percents 
above  lasing  threshold. 

The  gain  curve  in  Fig.  5  exhibits  two  peaks  and  a  dip  at 
fl  =  0.  For  T2<  7’ |  and  /; -  1  <€  1 ,  the  location  of  the  peaks  is 
given  by  (see  Appendix  A) 

~  ^  ^37 yr2'  ^ 

Note  that  Eq.  (23)  gives  a  smaller  frequency  than  Eq.  (19).  In 
addition,  the  splitting  in  Fig.  5  scales  like  the  square  root  of 
the  Rabi  frequency.  The  cw  solution  is  destabilized  when  the 
cavity  admits  a  mode  for  which  gill)  in  Fig.  5  is  positive. 

B.  Numerical  results 

The  moment  the  cw  solution  is  destabilized,  studying  Eqs. 
(7)-(10)  requires  a  numerical  simulation.  The  results  of  such 
a  simulation  are  shown  in  Fig.  6.  The  parameters  are  given  in 
Table  I,  with  the  only  difference  that  y=  0  was  used. 

Figure  6  shows  a  clear  pattern  in  the  spectrum.  This  pat¬ 
tern  appears  only  after  very  long  averaging  (order  of  a  mi¬ 
crosecond)  of  the  spectrum.  Such  averaging  is  appropriate. 


since  similar  averaging  occurs  when  the  spectra  are  mea¬ 
sured.  The  pattern  depends  on  T2  and  f0,  but  very  weakly 
depends  on  7).  A  similar  pattern  occurs  in  the  experimental 
spectra.  However,  we  were  not  able  to  trace  its  origin. 

The  combined  effect  of  SHB  and  a  saturable  absorber  is 
demonstrated  in  Fig.  7.  One  can  observe  two  effects.  First  is 
spectral  broadening  due  to  SHB,  similarly  to  Fig.  6.  Second, 
however,  is  the  appearance  of  a  splitting  in  the  spectrum. 
This  splitting  is  roughly  equal  to  twice  the  Rabi  frequency. 
Note,  however,  that  in  contrast  to  a  ring  laser,  in  a  Fabry- 
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FIG.  7.  Results  of  numerical  simulations  of  the  spectra  based  on 
the  Maxwell-Bloch  equations  including  a  saturable  absorber  and 
spatial  hole  burning  for  different  values  of  the  current  density  nor¬ 
malized  to  the  threshold  value. 
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FIG.  8.  (Color  online)  Scanning  electron  microscope  (SEM)  im¬ 
age  of  the  cross  section  of  a  10  fi m  wide  ridge  QC  laser. 

Perot  laser  “the  Rabi  frequency”  is  not  a  perfectly  well- 
defined  quantity.  Since  a  standing  wave  is  formed  in  the 
cavity,  the  electric  field  and  thus  the  Rabi  frequency  depend 
on  the  position  in  the  cavity.  This  dependence  is  even  stron¬ 
ger  when  the  mirrors  have  a  relatively  low  reflection  coeffi¬ 
cient,  since  the  field  amplitude  even  more  strongly  depends 
on  position. 

In  the  previous  section  we  saw  that  a  saturable  absorber 
lowers  the  threshold  of  the  RNGH  instability.  Here  we  see 
that  SHB  suppresses  the  central  peak  seen  in  Fig.  4  and 
replaces  it  by  a  minimum  in  the  spectrum.  In  the  next  sec¬ 
tions  we  see  that  Fig.  7  agrees  well  with  experimentally  mea¬ 
sured  spectra,  at  least  for  devices  where  the  RNGH  behavior 
was  dominant. 

VI.  EXPERIMENTAL  RESULTS:  RIDGE  LASERS 

We  first  study  the  multimode  regimes  in  standard  ridge 
QCLs,  in  which  the  sidewalls  of  the  laser  ridges  are  covered 
by  a  thick  layer  of  electrically  plated  gold  contact,  see  Fig.  8 
for  a  typical  laser  cross  section.  The  active  region  of  the 
samples  tested  is  based  on  a  three -quantum-well  design  emit¬ 
ting  at  a  wavelength  X.  —  8  /rm  [18].  The  wafer  was  grown 
by  metal  organic  vapor  phase  epitaxy  (MOVPE).  Figure  9(a) 
shows  the  voltage-current  ( V-I )  and  light-current  (L-I)  char¬ 
acteristics  of  a  10  /r m  wide  laser  operated  in  cw  at  200  K, 
and  Fig.  9(b)  shows  the  corresponding  optical  spectra.  The 
laser  was  cleaved  into  a  2  mm  long  bar  and  soldered  with 
indium  onto  a  copper  heat  sink.  The  optical  power  was  mea¬ 
sured  by  an  thermal  head  power  meter  with  a  collection  ef¬ 
ficiency  of  nearly  100%.  The  spectra  were  measured  by  a 
Nicolet  Fourier  transform  infrared  spectrometer  (FTIR) 
equipped  with  a  deuterated  triglycine  sulfate  (DTGS)  detec¬ 
tor. 

As  shown  in  Fig.  9(b),  the  laser  spectrum  is  single  mode 
close  to  laser  threshold.  It  broadens  and  splits  into  two  sepa¬ 
rated  humps  as  the  pumping  current  increases.  The  separa¬ 
tion  between  the  two  peaks  of  the  two  humps  increases  lin¬ 
early  with  the  square  root  of  the  collected  output  power  from 
one  facet,  as  shown  in  Fig.  9(c).  The  Rabi  angular  frequency 
can  be  calculated  from  the  collected  output  power,  using  the 
formula  HRabi  = =  /rV2n/ave/(ce)/^.,  where  /j,  is  the 
electron  charge  times  the  matrix  element  of  the  laser  transi¬ 
tion  (=1.9  nm  for  this  particular  device).  7ave  is  the  average 
intracavity  intensity  in  the  gain  region,  which  can  be  derived 
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from  the  measured  output  power  [19].  For  all  values  of  in¬ 
tensity  corresponding  to  the  spectra  reported  in  Fig.  9(b), 
nRabj/27r  was  calculated,  multiplied  by  a  factor  of  2  and  then 
added  to  Fig.  9(c)  (solid  line).  Reasonably  good  agreement  is 
found  between  the  experimental  splitting  and  twice  the  esti¬ 
mated  Rabi  frequency.  The  error  bars  of  the  spectrum  split¬ 
tings  come  from  the  uncertainty  in  determining  the  exact 
position  of  the  peaks,  which  is  the  full  width  at  half  maxi¬ 
mum  (FWHM)  of  the  humps.  As  mentioned  previously  in  the 
theoretical  section,  the  RNGH  instability  predicts  that  large 
intracavity  intensity  will  result  in  parametric  gain  at  frequen¬ 
cies  detuned  from  the  maximum  of  the  gain  curve  by  the 
Rabi  frequency.  The  measured  spectra  thus  show  clear  indi¬ 
cation  of  the  RNGH  instability  in  ridge  QCLs. 

The  lowering  of  the  RNGH  instability  threshold  in  our 
QCLs  is  due  to  the  presence  of  a  saturable  absorber.  This 
phenomenon  is  demonstrated  analytically  in  the  theoretical 
section.  Such  a  saturable  absorption  mechanism  in  our  ex¬ 
periments  is  likely  to  come  from  Kerr  lensing,  caused  by  a 
nonlinear  (i.e.,  intensity  dependent)  refractive  index  n2I  in 
the  active  region  [5].  As  the  light  intensity  increases,  the 
mode  becomes  more  confined  in  the  plane  transverse  to  the 
propagation  direction,  and  the  corresponding  net  modal  gain 
also  increases.  The  reason  is  twofold:  Lirst,  the  mode  over¬ 
laps  more  with  the  active  region,  leading  to  a  larger  modal 
gain  (this  mechanism  is  often  called  “soft  Kerr-lensing”). 
Second,  the  overlap  with  the  metal  contacts  is  reduced,  lead¬ 
ing  to  smaller  losses.  Thus  the  metal  acts  like  a  saturable 
absorber.  A  detailed  analysis  is  given  in  Appendix  C. 

The  same  RNGH  splitting  in  spectra  is  observed  in  many 
different  devices,  from  wafers  grown  by  both  molecular 
beam  epitaxy  (MBE)  and  MOVPE.  Figure  10  shows  the 
spectra  of  a  laser  fabricated  from  a  MBE-grown  wafer  with 
the  same  active  region  design  as  in  Fig.  9,  taken  in  continu¬ 
ous  wave  at  77  K.  The  spectrum  starts  from  single  mode 
close  to  the  laser  threshold,  and  the  mode  hops  at  a  pumping 
ratio  of  /Qth=  1.66,  where  there  is  a  corresponding  kink  in 
the  ( L-I )  curve  at  330  mA.  At  590  mA,  there  is  another  kink 
in  the  {L-I)  curve,  and  a  corresponding  abrupt  shift  in  the 
spectrum  associated  with  mode  hopping.  Mode  hopping  and 
its  relationship  to  kinks  in  the  {L-I)  curves  are  common  phe¬ 
nomena  in  semiconductor  lasers  [20-22],  but  its  cause  in  QC 
lasers  has  not  yet  been  studied  and  is  beyond  the  focus  of 
this  paper.  After  the  spectrum  broadens,  it  forms  two  sepa¬ 
rated  humps  whose  peaks  shift  apart  with  increased  optical 
power,  similar  to  the  MOVPE-grown  lasers. 

Characterization  of  the  pulse  emission  was  done  using  the 
method  of  second-harmonic  interferometric  autocorrelation. 
It  is  based  on  a  Michelson  interferometer  in  which  the  input 
beam  is  split  into  two  and  one  of  them  is  delayed  by  t.  In  the 
standard  setup,  once  recombined,  the  two  pulses  are  sent 
collinearly  first  into  a  nonlinear  crystal,  and  then  a  filter, 
which  allows  only  the  second-harmonic  generation  (SHG) 
component  to  be  detected.  One  can  test  if  it  originates  from  a 
stationary  isolated  pulse  train  from  the  ratio  between  the  in¬ 
terference  maximum  and  the  background  (see  Appendix  B). 
The  pulse  duration  can  also  be  determined.  However,  due  to 
the  extremely  low  SHG  conversion  efficiency  of  nonlinear 
crystals  in  the  mid-IR,  the  conventional  setup  is  not  feasible. 
To  overcome  this  problem,  we  used  a  two-photon  quantum 


053804-7 


GORDON  el  al. 


PHYSICAL  REVIEW  A  77,  053804  (2008) 


E, 

i 

CL 


CL 

o 


Current  (A) 


o 

Q_ 


Q. 

O 


_ ^^AAAAAAAAAAAAAAAA/iA^liWAAAAAAAAAAAA^VVLAAAAAAAAAAAAAA^^^^^^Ajv^^ _ 

1.88 

1.75 

1.63 

1.50 

1.38 

1.25 

1.13 

1.10 

1.08 

JL  ~ . 

j/jth=1.°5  1 

_ LJD _ I _ 1 _ I _ 1 _ I _ 1 - V-I _ 1 _ 1 _ 1 _ 1 _ 

(b) 


1280 


1290 


1300 


1310 


1320 


1330 


N 

I 

H 

& 

c 

<D 

=3 

W 

0) 


1340 


Wave  Number  (cm  ) 


Power 1/2  (mW1/2) 


FIG.  9.  (Color  online)  (a)  V-I  and  L-l  curves  and  (b)  optical  spectra  vs  pumping  ratio  (Jl  jth)  above  threshold  obtained  in  cw  at  200  K  of 
a  10  /j.m  wide  ridge  laser  emitting  at  8  gm  (wafer  No.  2721).  (c)  Spectral  splitting  and  twice  the  Rabi  frequency  HRabi/(2Tr)  vs  square  root 
of  output  power  collected  from  a  single  laser  facet.  The  different  quantities  reported  on  the  graph  were  deduced  from  the  experimental  data 
shown  in  (b).  The  dashed  line  is  a  least-squares  linear  fit  of  the  data. 


well  infrared  photodetector  (QWIP)  [23,24]  instead  of  using 
a  nonlinear  crystal  plus  a  linear  detector.  The  energy  diagram 
of  one  period  of  the  multiquantum  well  detector  under  bias  is 
shown  in  the  inset  of  Fig.  1 1 .  The  first  three  electronic  states 
are  nearly  equidistant  in  energy.  When  electrons  in  the  doped 
quantum  wells  absorb  two  photons  simultaneously  and  the 
detector  is  biased  (1-3  V),  a  photocurrent  is  generated  and 
the  signal  can  be  detected  by  use  of  a  preamplifier  and  a 
lock-in  amplifier.  This  experimental  setup  is  shown  in  Fig. 

11.  The  second-order  interferometric  autocorrelation  trace  of 
the  MBE-grown  device  mentioned  above  is  shown  in  Fig. 

12.  Interference  fringes  are  observed  when  the  delay  time 
from  one  arm  of  the  autocorrelator  is  equal  to  multiples  of 
the  cavity  round-trip  time.  The  ratio  between  the  maximum 
of  the  interference  fringes  and  the  background  is  smaller  than 
8,  and  the  autocorrelation  trace  has  some  features  between 
the  cavity  round-trip  times,  indicating  that  the  multimode 
regime  observed  in  this  device  is  not  stable  mode  locking 


with  a  single  pulse  per  round-trip.  The  meaning  of  the  auto¬ 
correlation  traces  is  discussed  in  Appendix  B,  and  simulated 
traces  are  given  therein. 

In  addition  to  second-order  autocorrelation,  the  micro- 
wave  spectrum  of  the  laser  output  was  also  measured  with  an 
ultrafast  QWIP  [25]  whose  bandwidth  is  52  GHz.  The  laser 
output  is  sent  directly  to  the  ultrafast  QWIP,  and  the  resulting 
photocurrent  is  displayed  in  a  spectrum  analyzer.  Figure  12 
shows  the  beat  note  signal  of  the  MOVPE-grown  sample  at 
pump  current  800  mA  at  77  K.  A  steady  peak  with  a  FWHM 
of  13  MHz  at  22.01  GHz,  which  corresponds  to  the  cavity 
round-trip  frequency  of  the  2  mm  long  laser  (background 
refractive  index  n  =  3.3),  is  observed  on  the  spectrum  ana¬ 
lyzer.  It  indicates  a  modulation  of  the  laser  output  at  the 
cavity  round-trip  frequency,  and  thus  at  least  partial  phase 
locking  between  the  longitudinal  modes:  The  phase  relation¬ 
ships  between  the  modes  are  stable  for  about  103  round-trips, 
as  inferred  from  the  peak  frequency-to-FWHM  ratio. 
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FIG.  10.  (Color  online)  (a)  V-I  and  L-l  curves  and  (b)  the  optical  spectra  vs  pumping  ratio  (// jth)  above  threshold  obtained  from  a  10  /rm 
wide  ridge  QCL  operated  in  cw  at  77  K  emitting  at  8  /j.m  (wafer  No.  2743).  Note  that  mode  hopping  occurs  while  there  is  a  jump  in  the  L-l 
curve. 


As  the  ridge  width  is  increased,  the  RNGH  instability  and 
the  Rabi  splitting  in  the  spectrum  are  suppressed.  Figure  13 
shows  the  spectra  of  a  device  processed  from  the  same  wafer 
as  the  one  in  Fig.  9,  with  the  difference  of  its  ridge  width 
being  increased  to  15  /tm.  The  spectra  do  not  broaden  much 
even  at  very  high  pumping  currents,  and  do  not  show  the 
Rabi  splitting.  As  discussed  in  Appendix  C,  the  effect  of  Kerr 
lensing  decreases  significantly  when  the  active  region  width 
is  increased.  Thus  this  is  strong  evidence  that  the  saturable 
absorption  needed  for  lowering  the  RNGH  threshold  is  pro¬ 
vided  by  the  Kerr-lensing  mechanism. 

VII.  EXPERIMENTAL  RESULTS:  BURIED 
HETEROSTRUCTURE  LASERS 

The  second  type  of  lasers  we  have  tested  are  the  so-called 
buried  heterostructure  lasers,  in  which  an  insulating  Fe- 
doped  InP  layer  is  regrown  after  etching  of  the  ridges.  A 


FIG.  1 1 .  (Color  online)  The  experimental  setup  of  a  two-photon 
autocorrelation  measurement.  (Inset)  Conduction  band  diagram  of 
the  two-photon  QWIP  showing  three  equidistant  energy  levels. 


thick  layer  of  electrically  plated  gold  is  deposited  on  top  of 
the  ridges  as  the  top  contact  layer  after  the  InP  regrowth. 
Figure  14  shows  the  cross  section  of  such  a  device. 

The  active  region  of  the  buried  heterostructure  lasers 
tested  is  based  on  a  four-quantum-well  design,  which  relies 
on  a  double  phonon  resonance  to  achieve  population  inver¬ 
sion  [26].  Figure  15(a)  shows  the  V-I  and  L-I  characteristics 
of  a  laser  with  an  active  region  width  of  3  /jl m  and  wave¬ 
length  8.38  /rm  operated  in  cw  at  room  temperature,  and 
Fig.  15(b)  shows  its  optical  spectra.  The  spectra  show  Rabi 
splitting  similar  to  the  ridge  laser  case,  indicating  the  RNGH 
instability  in  this  narrow  buried  heterostructure  laser.  The 
spectral  splitting  and  twice  the  Rabi  frequency  are  plotted 
against  the  square  root  of  the  collected  output  power  in  Fig. 
15(b).  Good  agreement  is  found  between  the  experimental 
splitting  and  twice  the  calculated  Rabi  frequency. 

The  second-order  interferometric  autocorrelation  trace  of 
the  device  when  it  is  pumped  five  times  above  threshold  at 
80  K  is  shown  in  Fig.  16.  The  ratio  between  the  maximum  of 
the  interference  fringes  and  the  background  is  close  to  8  to  3, 
similar  to  the  ridge  laser  case.  In  addition,  there  are  smaller 
interference  fringes  within  one  cavity  round-trip  period.  Both 
features  indicate  no  stable  pulsation  from  the  laser  (see  Ap¬ 
pendix  B).  The  microwave  spectrum  of  the  laser  output 
shows  a  steady  peak  with  FWHM  of  700  kHz  at  15.018 
GHz.  This  is  shown  in  the  inset  of  Fig.  16.  This  indicates  that 
the  phase  relationship  between  the  modes  are  maintained  for 
about  104  cavity  round-trips. 

As  in  the  ridge  laser  case,  the  lowering  of  the  threshold  of 
RNGH  instability  can  also  be  attributed  to  the  Kerr-lens  type 
saturable  absorption.  In  order  to  better  support  this  idea  ex¬ 
perimentally,  the  spectra  of  another  device  processed  from 
the  same  wafer  as  in  Fig.  15  but  with  a  wider  active  region 
(7.5  /rm)  are  measured.  Two-dimensional  waveguide  simu¬ 
lations  indicate  a  much  weaker  Kerr-lensing  effect  in  these 
QCLs,  due  to  the  much  larger  ratio  of  active  region  width  to 
wavelength  (see  Appendix  C).  The  measured  optical  spectra 
obtained  at  300  K  in  cw  mode  are  shown  in  Fig.  17(a).  The 
envelopes  of  the  spectra  consist  of  multiple  peaks  whose 
separation  is  independent  of  the  pumping  current.  The  spec- 
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FIG.  12.  (Color  online)  A 
second-order  interferometric  auto¬ 
correlation  trace  of  a  8  /jl m  wave¬ 
length  ridge  QC  laser  (wafer  No. 
2743)  under  the  condition  of 
RNGH  instability.  (Inset)  Micro- 
wave  spectrum  of  photocurrent 
generated  by  a  similar  laser  (wafer 
No.  2721)  under  the  condition  of 
RNGH  instability  (measured  with 
a  68  kHz  resolution  bandwidth). 


tral  signatures  qualitatively  agree  with  the  numerical  simula¬ 
tions  in  a  Fabry-Perot  cavity  without  a  saturable  absorber 
(7=0)  (Fig.  6). 

Further  increasing  the  active  region  width  tends  to  sup¬ 
press  both  the  instabilities  caused  by  RNGH  and  spatial  hole 
burning.  Figure  17(b)  shows  the  spectra  of  a  device  also 
processed  from  the  same  wafer  but  with  an  even  wider  active 
region  (10  /rm)  at  77  K.  The  spectra  do  not  broaden  much 
even  at  very  high  pumping  currents,  as  in  the  ridge  laser 
case.  Since  both  RNGH  instability  and  SHB  stem  from  non¬ 
linear  effects,  they  are  suppressed  when  the  intensity  of  the 
field  in  the  cavity  is  lower.  Moreover,  when  the  active  region 
width  increases,  higher  transverse  modes  are  excited.  Differ¬ 
ent  transverse  modes  have  different  propagation  constants  /3, 
and  thus  form  different  gain  gratings  which  tend  to  wash  out 
the  effect  of  spatial  hole  burning. 


FIG.  13.  Optical  spectra  vs  pumping  ratio  (Jl jth)  above  thresh¬ 
old  obtained  in  cw  at  77  K  with  a  15  /mi  wide  ridge  laser  emitting 
at  8  /zm  (wafer  No.  2721). 


VIII.  TEMPERATURE  EFFECTS 

Temperature  also  plays  an  important  role  in  the  shape  of 
the  spectrum.  To  illustrate  this  point,  we  now  present  the 
spectra  at  different  temperatures  for  the  \  =  8.38  /rm  buried 
heterostructure  laser  with  3  /tm  active  region  width  (the 
same  as  in  Fig.  15). 

From  the  spectra  in  Fig.  18,  it  is  clear  that  at  lower  tem¬ 
peratures  the  spectra  are  dominated  by  spatial  hole  burning, 
showing  multiple  peaks  independent  of  pumping  and  no  sig¬ 
nificant  Rabi  splitting.  As  the  temperature  increases,  the  Rabi 
splitting  becomes  more  evident  and  finally  the  RNGH  insta¬ 
bility  takes  over. 

The  effect  of  temperature  on  the  nature  of  the  multimode 
regime  in  QCLs  can  be  explained  in  part  by  carrier  diffusion, 
however,  this  effect  alone  seems  to  be  too  weak  to  provide  a 
complete  explanation.  The  lifetime  of  the  gain  grating  Tg  is 
defined  as  T~g  =T\X +Ak2D.  The  diffusion  coefficient  D  of  the 
gain  grating  is  proportional  to  the  temperature,  D  =  fxkBTI  q, 
where  /jl  here  is  the  carrier  mobility,  kB  is  the  Boltzmann’s 
constant,  T  is  the  temperature,  and  q  is  the  carrier  charge. 
Therefore  at  higher  temperature  carrier  diffusion  would  re- 


FIG.  14.  SEM  image  of  the  cross  section  of  a  buried  hetero¬ 
structure  QC  laser. 
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FIG.  15.  (Color  online)  (a)  V-I  and  L-I  curves  and  (b)  optical  spectra  vs  pumping  ratio  0/ 7th)  above  threshold  obtained  in  cw  at  300  K 
with  a  3  /x m  wide  ridge  laser  emitting  at  8.38  /rm  (wafer  No.  3251).  (c)  Spectral  splitting  and  twice  the  Rabi  frequency  HRabi/(2-7r)  vs 
square  root  of  output  power  collected  from  a  single  laser  facet.  The  dashed  line  is  a  least-squares  linear  fit  of  the  data. 


duce  spatial  hole  burning  and  thus  reveal  the  RNGH  insta¬ 
bility.  However,  with  an  upper  limit  for  the  mobility  of 
7000  (cm2/s)/V  and  A=2.25X104  cm-1  (which  corre¬ 
sponds  to  a  wavelength  in  vacuum  of  8.38  /rm),  4 k2D 
=  0.09  THz  at  77  K  and  4k2  D  =  0.4  THz  at  300  K,  both 
significantly  smaller  than  T^  —  O.h  THz  at  77  K  and  T\x 
~2  THz  at  300  K.  Thus  carrier  diffusion  is  unlikely  to  be 
the  reason  for  the  temperature  dependence.  Although  the 
temperature  effect  is  not  entirely  understood  to  us  at  this 
point,  one  possibility  is  temperature-dependent  saturable  ab¬ 
sorption.  The  QCL  injector  consists  of  many  energy  levels 
which  can  be  thermally  populated.  It  is  not  surprising  that  if 
any  two  higher  levels  in  the  injector  are  nearly  resonant  with 
the  laser  transition,  they  will  form  a  two-level  saturable  ab¬ 
sorber.  Therefore  in  this  case  saturable  absorption  is  stronger 
at  higher  temperature  and  makes  the  RNGH  instability  more 
easily  observable. 

IX.  CONCLUSION 

This  paper  provides  a  thorough  account  of  different  mul¬ 
timode  regimes  in  QCLs.  It  was  found  that  two  key  mecha¬ 


nisms  which  govern  the  multimode  regimes  in  QCLs  are  a 
coherent  instability  similar  to  the  RNGH  instability  and  spa¬ 
tial  hole  burning.  The  former  is  enhanced  due  to  the  large 
dipole  moment  /x  of  the  laser  transition,  which  results  in  a 
large  Rabi  frequency  compared  to  the  relaxation  rates.  Thus 
the  Rabi  splitting  can  be  resolved  by  the  comb  of  modes 
supported  by  the  cavity.  SHB  is  enhanced  because  carrier 
diffusion  is  slower  than  the  gain  recovery,  and  thus  leaves 
the  gain  grating  intact.  Note  that  in  conventional  semicon¬ 
ductor  lasers  the  RNGH  instability  is  not  observed  because 
typical  Rabi  frequencies  are  much  smaller  than  the  phase 
relaxation  rate  1IT2 .  SHB  in  diode  lasers  is  not  so  readily 
observable  because  diffusion  occurs  on  a  time  scale  compa¬ 
rable  to  the  recombination  time,  i.e.,  the  second  term  on  the 
right-hand  side  of  Eq.  (22)  is  of  the  same  order  or  greater 
than  the  first  term  1  /  7j  =  1 09  s_1 

Due  to  the  fast  gain  recovery,  conventional  mode  locking, 
with  one  pulse  per  round-trip,  is  suppressed.  In  order  to 
achieve  conventional  mode  locking  in  QCLs,  one  needs  to 
design  a  QCL  with  a  slower  gain  recovery,  such  that  T \  be¬ 
comes  longer  than  or  comparable  to  the  cavity  round-trip. 
Efforts  in  this  direction  are  currently  underway. 
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FIG.  16.  (Color  online)  A 
second-order  interferometric  auto¬ 
correlation  trace  of  a  8  /zm  wave¬ 
length  buried  heterostructure  QC 
laser  (wafer  No.  3251)  under  the 
condition  of  RNGH  instability. 
(Inset)  Microwave  spectrum  of  the 
photocurrent  generated  by  the 
same  laser  (measured  with  a  68 
kHz  resolution  bandwidth). 
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APPENDIX  A:  LINEAR  STABILITY  ANALYSIS 
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This  appendix  elaborates  on  the  linear  stability  analysis  of 

the  Maxwell-Bloch  equations  in  a  ring  cavity  [Eqs.  ^  = _ 2 _ 

(13)— (15)]  and  in  the  Fabry-Perot  cavity  [Eqs.  (7)-(10)].  Athr2 


1.  Ring  cavity  Equations  (13)— (15)  are  then  rewritten  as 

We  begin  with  a  linear  stability  analysis  of  Eqs.  (13)— (15), 
with  €(£’)  =  f?0- ~)\E\2.  In  order  to  keep  the  expressions  from  =  _  Qji-ir)-  -(f0-  y\E\2)E 

becoming  too  cluttered  we  define  c  z  2 
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FIG.  17.  Optical  spectra  vs 
pumping  ratio  (Jl  jth)  above 
threshold  obtained  in  cw  at  300  K 
with  (a)  a  7.5  pm  wide  and  (b)  a 
10  pm  buried  heterostructure  la¬ 
ser  emitting  at  8  pm  (wafer  No. 
3251). 
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FIG.  18.  Optical  spectra  vs  pumping  ratio  (/'/ jth)  above  threshold  obtained  in  cw  of  a  buried  heterostructure  laser  emitting  at  8.38  /im 
at  different  temperatures:  (a)  80  K,  (b)  150  K,  (c)  200  K,  and  (d)  300  K. 


1  A  ~  77 

d.r>=  — E-  — , 

2  T2 


,  r  Pt o  A  ye»~  ^ 

dA= - -  —  +  i(E  7i  —  c.c.) . 

Tx  T2  Tx 


(A3) 


(A4) 


Setting  the  left-hand  sides  of  Eqs.  (A2)-(A4)  to  zero  we  find 
that  they  admit  a  steady  state  solution  of  the  form  E=E,  rj 
=  rj,  A  =  A.  E,  r),  A  are  constants  in  time  and  space  satisfying 


(A5) 


Adding  perturbations  SE,  Srj,  SA  to  the  steady  state  so¬ 
lution,  and  linearizing  Eqs.  (A2)-(A4)  with  respect  to  the 
perturbations,  one  obtains  a  set  of  linear  equations: 


1  -  ~  8r>, 

dtSyj=-(ASER+SAE)--^, 
z  12 


-  _  5A 

d,SA  =  -  T2AESEr  -2ESrj,-  — , 


(A8) 


(A9) 


-d,SER  =  -  dz8ER  +  Sr,:- (f0 -  3y£2)^  (A10) 

c  2 


and 


y=-(t0-yE2)E, 


p  +  l  =  \l-^j(l+E2T1T2). 


(A6) 


(A7) 


E  was  assumed  real,  since  it  can  be  always  chosen  so  without 
loss  of  generality. 


1  -  ~  £>77d 

d,SyR  =  --ASEI--^, 
L  i2 


d,SEj  =-dz SE ,  -<?%-(€  o  -  yE2)  — .  (A  1 2) 


SE, 

2 


(All) 


The  two  sets  of  equations,  (A8)-(A10)  and  (All)  and  (A12), 
are  decoupled,  and  can  be  thus  studied  independently.  Since 
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both  sets  of  equations  are  translationally  invariant,  plane 
waves  are  their  eigenfunctions.  For  Eqs.  (A8)-(A10)  we 
therefore  choose  S7fI{z,t)  =  8rjj(t)e'kz,  and  similarly  for  8A 


and  SEr.  The  stability  of  the  cw  solution  is  thus  determined 
by  the  eigenvalues  of  the  matrix 

^  '  777(^0- y£2) 

2 


M  = 


1  2 


i.\ 
-E 
2 


c  cl  1  3 

-  -  -f0  +  -yE~  -  l k 
n  n  n 


\-2E 


2  2 
yE^-CoE 


0 


7 


(A  13) 


If  all  eigenvalues  have  a  negative  real  part,  the  cw  solution  is 
stable. 

The  eigenvalues  of  M  can  be  easily  found  numerically. 
However,  it  is  always  more  enlightening  to  have  an  analyti¬ 
cal  expression.  To  this  aim  we  observe  that  of  the  three  fre¬ 
quency  parameters  in  M  [771,  771,  and  ^(€0+  yE2)\,  the  lat¬ 
ter,  which  is  the  inverse  cavity  photon  lifetime,  is  often  the 
slowest  one  in  QCLs.  According  to  Table  I,  ^(7= 0.05  THz, 
7y 1  =  2  THz,  and  Z^1  =  15  THz.  We  thus  derive  expressions 
for  the  eigenvalues  which  are  correct  to  first  order  in  f0  and 

r 

For  € o=0  and  y=0,  the  eigenvalue  with  the  greatest  real 
part  is  \0(k)  =  -ickl  n.  Putting  k(k)  =  \Q(k)  +  kl(k)  into  the 
characteristic  polynomial  of  M  and  equating  the  parts  which 
are  first  order  in  f0,  y,  and  X,,  one  arrives  at 

€qc  (nr1  +  onr2-2(p-i) 
max  *  2n(nr1  +  0(nr2  +  i)-(p-i) 

|  (p  -  l)yc  (nr,  +  Q(3nr2  +  2 f)  -  4 (P  - 1) 

2nTlT2  (ClTl  +  i){£lT2  +  i)-{p-\)  ’ 

(A  14) 


dAo  ~ 


A„  -  A0 


+  i(E*+r/+  +  E*_r/_-c.c.),  (A18) 


d,A2=  ±  i{E*+  r/_  —  rj+E_)  —  . 


(A19) 


Splitting  the  variables  with  the  tilde  into  their  real  and  imagi¬ 
nary  parts,  Eqs.  (A16)-(A19)  can  be  rewritten  as 

dM  +  #M-  =  |(A0  +  A  2)(SEr_  +  SE« )  +  E(8A0  +  SAR ) 


8rf+  +  8rj 


(A20) 


-{d,8ER+  +  dM)  =  -  dz(8ER  -  8Er)  +  (Srf++  Srf_) 
c 

(,  ~ 

-  —(8Er+SEr_),  (A21) 

d,8A0  =  2irj(8ER  +  8ER)  -  2 E(St?+  +  8rf_)  -  — , 

T\ 

(A22) 


d,8AR  =  ir,{8ER+  +  SEr)  -  E{8r}+  +  8rf_)  -  , 


(A23) 


dM  -  dM  =  -  2^0  -  Ao) ( 8E?+  -  8E'_)  +  8A2E 


8t)r  -  8rjR 


(A24) 


where  p= Ap/  Ath- 1  and  il  =  kcl n.  Taking  the  real  part  of  Eq. 
(A14)  one  obtains  Eq.  (17). 

The  position  of  the  maximum  of  the  gain  curve  is  inde¬ 
pendent  of  y,  and  to  first  order  in  T2!T\  is  given  by  Eq.  (19). 
The  gain  at  that  frequency  is  given  by 

g((^max))  =  ——[3 ip  -  1)  -  2^2p(]}  -  1)]  +  - {'  - . 

2 n  Ti  T{T2 

(A  15) 

The  threshold  for  the  instability  is  found,  to  leading  order  in 
77/77,  by  equating  Eq.  (A15)  to  zero.  This  yields  Eq.  (21). 

2.  Fabry-Perot  cavity 

Employing  the  same  transformation  as  in  Eq.  (Al),  Eqs. 
(7)-(10)  take  the  form 

n  ~  —  1  ~ 

—d,E±  =  -+-  3JE±  -  ir)±  -  ~(qE±,  (A16) 

c  2 

dtV±=l-(AaE±  +  AM)-r^,  (A17) 

Z  l2 


-( dM  -  dM)  =  -  dz(SE'+  +  8E‘_)  -M+-  M) 

c 

-^(8E,+  -8EI+),  (A25) 

d,8A[  =  ~iy( SE\  -  8E'_)  -  E( SpR  -  SpR_)  ~  • 

1g 

(A26) 

We  now  define  the  generalized  parity  operator  P,  which  re¬ 
flects  the  space  around  the  center  of  the  cavity  and  swaps  all 
+  and  -  quantities: 

PX±(z,t)=XT{L-z,t), 

where  X  is  any  of  the  quantities  rj,  E,  or  A,  with  an  R  or  7 
superscript.  Equations  (A20)-(A26)  commute  with  P,  and 
therefore  their  eigenmodes  can  be  chosen  to  be  eigenstates  of 
P  as  well.  Since  P1  is  the  unity  operator,  P  can  have  only 
“even”  or  “odd”  eigenstates,  with  +1  and  -1  eigenvalues, 
respectively. 
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We  now  observe  that  Eqs.  (A20)-(A23)  are  decoupled 
from  Eqs.  (A24)-(A26).  For  Eqs.  (A20)-(A23),  odd  solu¬ 
tions  are  trivial,  whereas  for  Eqs.  (A24)-(A26),  even  solu¬ 
tions  are  trivial.  Assuming  an  even  solution  for  Eqs. 
(A20)-(A23)  and  an  odd  solution  to  Eqs.  (A24)-(A26),  Eqs. 
(A20)-(A23)  are  reduced  to 


,  An  +  A t  E  ~  ~  Srf, 

dM~  2  ~SEr++  2(S\)+SAr)-  r'+, 

{All) 

-d,SER+  =  -  d,SER+  +  Srf+  —  —  SEr, 
c  ‘  2 

(A28) 

dtSA0  =  4;  rjSER  -  AE 8rj+  -  , 

7  i 

(A29) 

~n  ,  SA^ 

dtSA2  =  2i r)SER  -  2E Srj+  -  , 

1  8 

(A30) 

„  A0  —  A2  1  -  SrjR 

3, St—  2  2  dE+  +  2  SA2E  —  » 

(A31) 

-dtSEr+  =  -  dzSE'+  -  Srfl  -  SE{, 
c  2 

(A32) 

d,SA!2  =  -  2  i  r,SE\  -2ESt~—- 

7  e 

(A33) 

Similarly  to  the  discussion  around  Eq.  (A13),  the  stability 
of  Eqs.  (A27)-(A33)  is  studied  by  finding  the  eigenvalues  of 
the  matrix 
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for  Eqs.  (A31)-(A33).  E  is  now  related  to  p  via 


(A34) 


(A35) 


2TlT2  +  T2Tg 

The  matrix  (A34)  is  related  to  the  RNGH  instability,  whereas 
matrix  (A35)  is  related  to  spatial  hole  burning.  In  the  limit  of 


Tg—* 0,  as  well  as  in  the  case  Tg=Th  g((l)  is  identical  to  the 
expression  (17)  without  a  saturable  absorber  (y=0): 

{VLTx  +  i)VLT2-2(p-l) 

Wi+i)m2+i)-ip- 1). ' 

(A36) 

The  RNGH  instability  threshold  is  therefore  again  around 
p= 9.  Note,  however,  that  the  cw  solution  destabilizes  for 
much  smaller  p  due  to  spatial  hole  burning.  Indeed,  g(il) 
obtained  from  Eq.  (A35)  is  given  by 

3/(/  +  ar1)-(P-i) 

3(i  +  nr2)(i  +  nr2)-(p-i)_ 

(A37) 


g(H)  =  -  ~“Re 
2  n 


g(H)  =  -  ~"Re 
2  n 


for  Tg=Tx.  The  peak  of  the  gain  curve  is  obtained  at 


n 


2 

max 


1  Ip~1 

t1  \  3r,r2 


for  T2<  7)  and  p—  1  <g  1.  For  Tg= 0  spatial  hole  burning  does 
not  exists,  and  g(17)  obtained  from  Eq.  (A35)  is  given  by 


g(^): 


Ipc 
2  n 


Re 


2  i  +  ClT2 
i  +  ITT') 


(A38) 


which  is  never  positive. 


APPENDIX  B:  INTERFEROMETRIC  AUTOCORRELATION 


If  we  write  the  electric  field  as  a  function  of  time  as 
E(t)eu"’,  the  two-photon  interferometric  autocorrelation  is 
given  by 

r<x> 

7(r)=  |  E(!+T)eiu>T+E{t)\4dt 

J  —00 

r<x> 

=  [|£(f  +  r)|4  +  \E{t)\A  +  4|£(r  +  T)E(t)\2]dt 

J  —oo 

/•oo 

dt\E{t)\2E*{t)E{t  +  r)eiaT+  c.c. 

J  —oo 

dt\E(t+  T)\2E*(t)E(t+  T)eiaT+c.  c. 


+  2 


f 

{•CG 

J  — oc 


dt[E{t  +  T)E*(t)e2iu>r\2  +  c.c. 

) 

We  define  the  background  as 

h 


rco 

,  =  2  \E(t)\4dt. 

J  —00 


(Bl) 


(B2) 


From  Eq.  (Bl)  one  can  see  that  7(0)  =  8/fc.  We  now  assume 
that  E{t)  is  an  isolated  pulse,  and  for  simplicity’s  sake  we 
assume  that  E{t)  is  nonzero  only  over  an  interval  T.  Then  for 
t>T  all  terms  in  Eq.  (Bl)  which  include  both  t  and  t+r 
vanish.  Therefore  /(r>7’)  =  //).  In  other  words,  for  a  pulse 
one  has 
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FIG.  19.  (Color  online)  Simulated  interferometric  autocorrela¬ 
tion,  corresponding  to  the  same  parameters  as  in  Fig.  7.  The  value 
of  p  is  denoted  at  each  plot. 


m 


I(t>T) 


(B3) 


Now  let  us  assume  that  E{t)  is  a  complex  stationary  random 
process  whose  phase  at  each  point  in  time  is  uniformly  dis¬ 
tributed  over  the  interval  [0,27t].  Let  us  assume  for  simplic¬ 
ity  that  there  is  T  such  that  for  t>T  E(t )  and  E(t+  r)  are 
statistically  independent.  One  can  see  that  the  mean  value  of 
the  third,  the  fourth,  and  the  fifth  line  of  Eq.  (Bl)  vanish.  The 
second  line  of  Eq.  (Bl)  equals  3 Ib.  We  therefore  obtain  that 
for  a  random  process 


W)  _  8 

I(t>  T)  3' 


(B4) 


The  ratio  between  the  peak  and  the  background  in  the  auto¬ 
correlation  function  is  therefore  8:1  for  an  isolated  pulse,  that 
is  when  the  phases  of  all  modes  are  all  zero,  and  8:3  for 
modes  with  completely  random  phases. 


Figure  19  shows  simulated  interferometric  autocorrelation 
traces  for  the  same  parameters  as  in  Fig.  7.  For  p=  1  the 
phases  of  the  modes  are  random,  and  the  peak-to-background 
ratio  is  8:3.  When  the  saturable  absorber  is  absent,  y=0,  the 
ratio  in  the  autocorrelation  is  also  8:3.  When  the  pumping  is 
higher  (Fig.  19  below),  the  saturable  absorber  induces  some 
phase  relationships  between  the  modes,  and  the  peak-to- 
background  ratio  is  about  8:2.  Note  that  in  the  autocorrela¬ 
tion  in  Fig.  12,  the  ratio  is  slightly  greater  than  8:3.  The 
autocorrelation  also  has  a  nontrivial  structure,  with  a  peak  at 
half  the  cavity  round-trip  time.  We  were  not  able  to  repro¬ 
duce  this  structure  in  the  simulations. 


APPENDIX  C:  KERR  LENSING  IN  QCLs 

This  section  elaborates  on  the  analysis  of  the  Kerr  lensing 
effect  (including  soft  Kerr  lensing)  in  QCLs.  The  mechanism 
of  Kerr  lensing  can  be  understood  as  follows:  above  lasing 
threshold,  there  is  net  gain  in  the  active  region  which  com¬ 
pensates  the  mirror  losses.  Assuming  a  nonlinear  refractive 
index  A nNL=n2I  in  the  active  region,  it  causes  self-focusing 
of  the  transverse  mode,  which  results  in  a  net  increase  in  the 
modal  gain  due  to  an  increased  overlap  with  the  active  re¬ 
gion  and  a  decreased  overlap  with  lossy  waveguide  cladding. 
A  stronger  intensity  leads  to  stronger  Kerr  lensing,  thus 
forming  an  intensity-dependent  saturable  absorber. 

The  active  region  width  is  expected  to  play  an  important 
role  in  the  Kerr  lensing  effect.  The  narrower  the  active  re¬ 
gion,  the  less  confined  is  the  transverse  mode,  which  results 
in  a  bigger  increase  of  the  modal  gain  due  to  Kerr  lensing. 
Therefore  the  saturable  absorber  coefficient  y  is  expected  to 
be  larger  in  lasers  with  a  narrower  active  region. 

We  performed  FDTD  waveguide  simulations  using  the 
commercial  software  BeamPROP.  We  simulated  buried  het¬ 
erostructure  QCLs  (wafer  No.  3251)  with  various  active  re¬ 
gion  widths.  The  software  allows  us  to  assign  complex  re¬ 
fractive  indices  to  each  layer,  and  the  effective  modal  index 


FIG.  20.  (Color  online)  Simu¬ 
lation  results  of  the  modal  gain  vs 
the  nonlinear  refractive  index  n2I 
for  a  3  pm  wide  and  a  5  pm 
wide  active  region  in  wafer  No. 
3251. 
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is  calculated  by  the  software.  The  detailed  steps  of  our  Kerr- 
lensing  simulation  are  as  follows. 

(1)  We  first  assign  a  net  gain  (i.e.,  a  negative  imaginary 
part  of  nactiVe)  in  the  active  region,  making  the  total  modal 
gain  equal  to  the  total  mirror  loss.  In  our  case,  given  the 
reflectivity  r=0.21  at  our  laser  facets,  a  3  mm  long  laser  has 
a  total  mirror  loss  of  4  cm-1. 

(2)  We  assign  different  values  of  A nNL=n2I  to  the  real 
part  of  n  actiVe.  After  running  the  simulation,  we  get  the  net 
modal  gain  for  each  A nNL. 

(3)  We  plot  the  modal  gain  versus  A nNL.  For  example. 
Fig.  20  shows  the  plots  for  3  /rm  wide  and  5  ptm  wide 
active  regions. 

(4)  In  order  to  relate  to  our  theoretical  model,  the  change 
in  modal  gain  is  attributed  entirely  to  change  in  modal  losses. 

(5)  The  change  in  modal  loss  equals  y\E\1 2,  which  is  also 
proportional  to  the  intensity  I.  Thus  from  the  slope  of  the 
plot  of  modal  gain  vs  A nNL=n2I,  we  can  extract  the  ratio  of 

yin2- 

(6)  For  a  given  active  region  width,  we  choose  the  value 
of  y  that  brings  the  threshold  of  RNGH  instability  down 
close  to  our  experimental  data.  Then  we  can  determine  the 
value  of  n2  from  the  ratio  of  yln2. 

(7)  With  the  fixed  value  of  n2,  we  repeat  the  waveguide 
simulation  with  different  active  region  widths  and  plot  y  vs 
the  active  region  width. 

For  an  active  region  width  of  3  /xm,  the  RNGH  threshold 
is  dramatically  lowered  from  nine  times  the  lasing  threshold 
to  less  than  twice  the  lasing  threshold  (see  Fig.  15).  From 
Fig.  2,  the  required  y  is  about  2  X  10-9  cm/ V2,  from  which 
we  obtain  that  n2  is  about  2  X  10-8  cm2 4 5 6 7 8/ W.  The  estimated  y 
vs  active  region  width  is  plotted  in  Fig.  21.  As  the  active 
region  width  increases,  y  decreases,  and  when  the  active 
region  is  as  wide  as  7.5  /xm,  y  decreases  to  1 
X  1  (r9  cm/V2.  This  pushes  the  RNGH  instability  threshold 
to  above  twice  the  lasing  threshold.  This  may  explain  why 


FIG.  21.  (Color  online)  Simulation  results  of  the  saturable  ab¬ 
sorber  coefficient  y  vs  the  active  region  width  in  wafer  No.  3251 
assuming  n2  =  2  X  10-8  cm2/W. 

we  see  only  instabilities  resulting  from  spatial  hole  burning 
in  wider  buried  heterostructure  lasers  (see  Fig.  17). 

For  ridge  lasers,  the  lossy  gold  contact  on  the  sidewalls  of 
the  laser  ridges  will  make  the  saturable  absorber  even  stron¬ 
ger.  This  may  explain  the  absence  of  spatial  hole  burning  in 
ridge  lasers. 

There  are  several  possibilities  which  might  explain  the 
origin  of  n2  in  the  active  region.  The  complex  states  in  the 
injector  can  all  contribute  to  both  real  and  imaginary  parts  of 
y  21  at  the  laser  wavelength,  which  results  in  both  n2  and  y. 
The  laser  transition  itself  can  contribute  an  n2  as  big  as 
10-9  cm2/ W  [5].  y  should  be  enhanced  with  increasing  tem¬ 
perature  since  electrons  at  higher  temperatures  populate  a 
larger  number  of  excited  states,  from  which  absorption  can 
take  place;  and  although  not  very  significant,  there  is  always 
bulk  n2  of  the  material. 
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We  report  the  observation  of  a  coherent  multimode  instability  in  quantum  cascade  lasers  (QCLs),  which  is 
driven  by  the  same  fundamental  mechanism  of  Rabi  oscillations  as  the  elusive  Risken-Nummedal-Graham- 
Haken  (RNGH)  instability  predicted  40  years  ago  for  ring  lasers.  The  threshold  of  the  observed  instability  is 
significantly  lower  than  in  the  original  RNGH  instability,  which  we  attribute  to  saturable-absorption  nonlin¬ 
earity  in  the  laser.  Coherent  effects,  which  cannot  be  reproduced  by  standard  laser  rate  equations,  can  play 
therefore  a  key  role  in  the  multimode  dynamics  of  QCLs,  and  in  lasers  with  fast  gain  recovery  in  general. 
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The  fundamental  coherent  mechanism  that  can  destabilize 
a  single-mode  laser  was  predicted  in  the  early  1960s  [1]  and 
was  later  extended  to  multimode  lasers  [2,3]  where  it  became 
known  as  the  Risken-Nummedal-Graham-Haken  (RNGH) 
instability.  These  instabilities  became  classic  landmarks  for 
the  general  field  of  nonlinear  dynamics  [4,5]  because  they 
emerge  in  conceptually  the  simplest  laser  model,  which  in 
the  single-mode  case  was  shown  to  be  equivalent  to  the  Lor- 
entz  model  of  deterministic  chaos  [6].  Another  feature  that 
makes  these  instabilities  so  interesting  and  unique  is  their 
coherent  nature  that  involves  the  polarization  of  the  medium 
as  a  nontrivial  dynamical  variable.  Most  other  physical 
mechanisms  that  can  drive  a  laser  from  a  single-mode  to  a 
multimode  regime,  such  as  spatial  and  spectral  hole  burning, 
0  switching,  and  saturable  absorption  [7,8],  can  be  ad¬ 
equately  described  within  the  standard  rate  equation  formal¬ 
ism,  in  which  the  polarization  of  the  active  medium  is  adia- 
batically  eliminated.  Both  the  single  mode  [1]  and  the 
multimode  [2,3]  instabilities  cannot  be  explained  by  the  rate 
equations.  Such  coherent  effects  can  be  only  observed  when 
the  polarization  is  driven  faster  than  or  on  a  time  scale  com¬ 
parable  to  the  dephasing  time  T  2  [9]. 

The  origin  of  the  two  coherent  phenomena  mentioned 
above  is  the  oscillation  of  the  population  inversion  at  the 
Rabi  frequency  QRabi  that  takes  place  when  the  intracavity 
laser  intensity  becomes  large.  This  results  in  a  modification 
of  the  gain  spectrum  and  the  emergence  of  sidebands  sepa¬ 
rated  from  the  maximum  of  the  gain  curve  by  an  amount 
corresponding  to  the  Rabi  frequency.  These  sidebands  can  be 
regarded  as  a  manifestation  of  parametric  gain.  The  instabil¬ 
ity  sets  in  when  the  intracavity  power  is  sufficiently  large: 
the  Rabi  angular  frequency  ilRabi  has  to  be  greater  than  the 
relaxation  time  scales  of  the  gain  medium  [more  precisely. 


*Electronic  address:  kaertner@mit.edu 
'Electronic  address:  capasso@deas.harvard.edu 


nRabi  is  sufficiently  greater  than  (T1T2)~1/2,  where  T  t  is  the 
gain  relaxation  time].  The  instability  threshold  is  often  called 
the  second  laser  threshold  due  to  its  universal  nature. 

Pioneering  theoretical  works  stimulated  extensive  experi¬ 
mental  studies  that  finally  resulted  in  the  observation  of  the 
Lorenz-type  chaos  in  a  far-infrared  single-mode  laser  [10]. 
However,  despite  almost  40  years  of  efforts,  the  experimen¬ 
tal  demonstration  of  the  multimode  RNGH  instability  has 
remained  somewhat  controversial  [11-16]. 

In  lasers  with  long  gain  recovery  compared  to  the  cavity 
round-trip  time,  the  instability  caused  by  a  saturable  absorber 
can  often  lead  to  mode  locking  [8].  When  the  gain  recovery 
time  is  short  compared  with  the  cavity  round-trip  time,  it  is 
usually  assumed  that  the  laser  dynamics  becomes  very  primi¬ 
tive  and  uninteresting  (so-called  class  A  laser).  In  this  case 
mode  locking  is  impossible  according  to  conventional  theory, 
and  the  relaxation  oscillation  frequency  becomes  purely 
imaginary  [17].  Surprisingly,  as  we  show  in  this  Rapid  Com¬ 
munication,  it  is  under  these  conditions  that  the  elusive 
RNGH  instability  can  be  observed.  We  show  that  quantum 
cascade  lasers  (QCLs)  are  uniquely  suited  for  studying  these 
coherent  effects  which,  along  with  spatial  hole  burning 
(SHB),  become  a  key  factor  in  dictating  the  dynamics  of  the 
laser. 

QCLs,  because  they  are  based  on  ultrafast  tunneling  and 
phonon-limited  intersubband  transitions,  belong  to  the  class 
of  lasers  which  have  an  extremely  fast  gain  recovery,  in  the 
range  of  a  few  ps  [18].  Recent  experiments  showed  indeed 
that  the  gain  recovers  within  a  few  ps,  which  is  approxi¬ 
mately  an  order  of  magnitude  shorter  than  the  cavity  round- 
trip  time  [19].  Since  its  invention  in  1994,  QCLs  have  un¬ 
dergone  tremendous  improvement  [20].  Recent  development 
of  low  loss,  high  power  QCLs  [21,22]  enables  the  study  of 
previously  under  investigated  aspects,  such  as  the  richness  of 
the  optical  spectrum  and  the  ultrafast  dynamics  of  these  de¬ 
vices.  In  Ref.  [23],  strong  evidence  of  self-pulsations  at  the 
cavity  round-trip  frequency  was  reported  in  QCLs,  in  par- 


1050-2947/2007/75(3)7031802(4) 


031802-1 


©2007  The  American  Physical  Society 


RAPID  COMMUNICATIONS 


WANG  etal. 


PHYSICAL  REVIEW  A  75,  031802(R)  (2007) 


ticular  a  large  broadening  of  the  spectrum  above  the  thresh¬ 
old  of  this  instability  was  observed.  However,  no  detailed 
pulse  characterization  was  provided.  The  technological  po¬ 
tential  of  QCLs  calls  for  a  better  understanding  of  the  inter¬ 
play  of  various  instabilities  in  the  parameter  regime  dictated 
by  these  lasers.  Moreover,  the  Rabi  frequency  in  QCLs  at  the 
power  levels  of  a  few  hundred  mW  is  of  the  order  of  a  few 
THz,  much  larger  than  the  spacing  of  Fabry-Perot  modes. 
Therefore  coherent  effects  should  be  easily  observable  in 
QCLs. 

In  this  Rapid  Communication  we  present  a  clear  experi¬ 
mental  demonstration  of  a  coherent  instability,  driven  by  the 
same  mechanism  as  the  RNGH  instability.  It  is  identified  in 
the  most  direct  manner,  by  demonstrating  in  the  optical  spec¬ 
trum  of  QCLs  a  splitting  corresponding  to  twice  the  Rabi 
frequency. 

The  instability  observed  differs  in  some  aspects  from  the 
original  RNGH  instability  [2,3].  The  threshold  of  instability 
can  be  as  low  as  a  few  tens  of  percent  above  the  laser  thresh¬ 
old,  as  shown  in  Fig.  1(a).  In  addition,  the  pure  RNGH  in¬ 
stability  typically  gives  rise  to  spectra  with  one  central  mode 
and  two  sidebands  separated  from  it  by  the  Rabi  frequency, 
whereas  in  our  experiments  we  observed  two  peaks  only, 
similarly  to  Ref.  [11].  However,  the  mechanism  of  the  insta¬ 
bility  is  the  same  in  essence,  namely  the  Rabi  oscillations  of 
the  population  inversion  due  to  coherent  effects.  The  differ¬ 
ences  from  the  RNGH  instability  as  it  occurs  in  ideal  condi¬ 
tions  [2,3]  can  be  attributed  to  the  presence  of  saturable  ab¬ 
sorption  and  SHB. 

The  QCLs  studied  were  fabricated  from  two  different  wa¬ 
fers  (wafer  nos.  3251  and  3252)  grown  by  metalorganic  va¬ 
por  phase  epitaxy.  The  devices  were  processed  into  buried 
heterostructure  lasers,  in  which  an  insulating  Fe-doped  InP 
layer  is  regrown  after  etching  of  the  ridges  [21,22].  The  ac¬ 
tive  region  of  all  the  samples  tested  is  based  on  a  four- 
quantum-well  design,  which  relies  on  a  double  phonon  reso¬ 
nance  to  achieve  population  inversion  [24].  Note,  however, 
that  the  multimode  operation  described  in  the  present  paper 
was  also  observed  with  lasers  based  on  so-called  three- 
quantum-well  designs  [18].  Figure  1(a)  shows  the  optical 
spectra  of  a  laser  operated  in  continuous  wave  (cw)  at  room 
temperature.  The  active  region  of  this  laser  is  3  /um  wide  and 
its  emission  wavelength  is  close  to  8.47  /xm  (wafer  no. 
3251).  The  laser  was  cleaved  into  a  3-mm-long  bar  and  sol¬ 
dered  with  indium  onto  a  copper  heat  sink.  The  spectra  were 
measured  by  a  Fourier  transform  infrared  spectrometer 
(FTIR)  equipped  with  a  deuterated  triglycine  sulphate 
(DTGS)  detector. 

As  shown  in  Fig.  1(a),  the  laser  spectrum  is  single  mode 
close  to  threshold  and  broadens  as  the  pumping  current  in¬ 
creases,  splitting  into  two  separated  humps.  The  difference 
between  the  weighted  centers  of  the  two  peaks  increases  lin¬ 
early  as  a  function  of  the  square  root  of  the  collected  output 
power  from  one  facet,  as  shown  in  Fig.  1(b)  (square  dots 
with  the  dashed  line  as  its  best  fit).  The  Rabi  angular  fre¬ 
quency  QRabi  can  be  easily  calculated  using  the  formula 
flRabi=/u£ th  =  jLLpnl avc/ (c e) / h,  where  /i  is  the  electron 
charge  times  the  matrix  element  of  the  laser  transition 
(=2.54  nm).  /ave  is  the  average  intracavity  intensity  in  the 
gain  region,  which  can  be  derived  from  the  measured  output 
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FIG.  1.  (Color  online)  (a)  Optical  spectra  vs  pumping  ratio 
O'/Jth)  above  threshold  obtained  in  cw  at  300  K  with  a  3-/um-wide 
buried  heterostructure  lasers  emitting  at  8.47  /tm.  For  1  <7  /yth 
<1.2  the  spectra  are  identical  to  j  / )th=  1 .2.  (b)  Spectral  splitting 
and  twice  the  Rabi  frequency  flRabi/(2Tr)  vs  square  root  of  output 
power  collected  from  a  single  laser  facet.  The  different  quantities 
reported  on  the  graph  were  deduced  from  the  experimental  data 
shown  in  (a).  The  dashed  line  is  a  least-square  linear  fit  of  the  data. 

power.  C  is  the  speed  of  light  in  vacuum  and  n  is  the  back¬ 
ground  refractive  index.  For  all  the  values  of  the  intensity 
corresponding  to  the  spectra  reported  in  Fig.  1(a),  llRahl  was 
calculated,  multiplied  by  a  factor  2  and  then  added  to  Fig. 
2(b)  (solid  line).  A  very  good  agreement  is  found  between 
the  experimental  splitting  and  twice  the  estimated  Rabi  fre¬ 
quency.  Both  curves  fall  indeed  well  within  the  error  bars 
[25].  As  mentioned  before,  the  theory  behind  the  RNGH  in¬ 
stability  predicts  that  the  large  intracavity  intensity  will  result 
in  parametric  gain  at  frequencies  detuned  from  the  maximum 
of  the  gain  curve  by  the  Rabi  frequency.  The  agreement  men¬ 
tioned  above  is  thus  a  strong  indication  of  the  RNGH  insta¬ 
bility  in  QCLs. 

In  order  to  better  understand  the  experimental  spectra  of 
the  QCLs  presented  in  Fig.  1(a),  we  use  a  simple  model 
based  on  the  standard  one-dimensional  Maxwell-Bloch  equa¬ 
tions  [9],  where  the  active  medium  is  described  by  an  “open” 
two-level  system  [26].  However,  contrary  to  the  standard 
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FIG.  2.  Results  of  numerical  simulations  of  the  spectra  based  on 
the  Maxwell-Bloch  equations  including  a  saturable  absorber  and 
spatial  hole  burning  for  different  values  of  the  current  density  nor¬ 
malized  to  the  threshold  value. 


unidirectional  Maxwell-Bloch  equations,  we  allow  the  elec¬ 
tromagnetic  field  to  propagate  in  both  directions.  The  waves 
traveling  in  the  two  directions  are  coupled,  as  they  share  the 
same  gain  medium.  This  gives  rise  to  SHB  [7]:  The  standing 
wave  formed  by  a  cavity  mode  imprints  a  grating  in  the  gain 
medium  through  gain  saturation.  As  a  result,  other  modes 
become  more  favorable  for  lasing,  and  a  multimode  opera¬ 
tion  is  triggered. 

In  the  slowly  varying  envelope  approximation,  the  equa¬ 
tions  read 

n  _  kl\lur  1 

~dtE±=  +  dzE±-l^—Irf±-  -£(£+,£_)£+,  (la) 

dtV±=i^j(\)E±  +  ^tE+)-^’  (lb) 

2  n  T  2 

dtA0=-^~ - +  -JL(E  +  7J+  +  E_7J_-C.C.),  (lc) 

/  1  ft 

dAt=  ±iT(E*+V-~  V+E-)~y-  (ld) 

ft  I  i 

The  +  and  —  subscripts  label  the  two  directions  of  propaga¬ 
tion.  £  and  77  are  the  slowly  varying  envelopes  of  the  electric 
field  and  the  polarization,  respectively.  The  actual  electric 
field  and  polarization  are  obtained  by  multiplying  E  and  77  by 
eltot  (a>  is  the  optical  resonance  frequency)  and  taking  the  real 
part.  The  position-dependent  inversion  is  written  as  the  sum 
of  the  three  terms,  A0,  A+e2ikz,  and  A^e~2ikz,  where  (A})* 
=  A2.  The  inversion  is  thereby  represented  by  two  slowly 
varying  functions  (A0  and  AJ),  and  e±2lkz  gives  the  fast  varia¬ 
tion  in  space.  All  the  quantities  mentioned  so  far  are  func¬ 
tions  of  space  Z  and  time  t. 


£(£  +  ,£_)  is  the  loss  in  the  cavity  (not  including  the  mirror 
loss),  which  is  allowed  to  be  nonlinear  and  dependent  on  the 
intensity.  In  this  work  we  assume 


£(£+,£_)  =  €0-r(|£+|2+|£_|2),  (2) 

where  £0  is  the  linear  loss  and  y  is  the  self-amplitude  modu¬ 
lation  coefficient  characterizing  the  nonlinear  (saturable)  part 
of  the  loss.  Such  a  saturable  absorption  mechanism  can  come 
from  Kerr  lensing  [8,23],  caused  by  a  nonlinear  refractive 
index  n2  in  the  active  region.  As  the  intensity  increases,  the 
mode  is  more  confined  in  the  plane  transverse  to  the  propa¬ 
gation  direction,  and  the  net  gain  it  experiences  is  greater. 
The  reason  is  twofold:  First,  the  mode  overlaps  more  with 
the  active  region,  leading  to  a  larger  modal  gain  (this  mecha¬ 
nism  is  often  called  “soft  Kerr  lensing”  [27]).  Second,  the 
overlap  with  the  metal  contacts  is  smaller,  leading  to  smaller 
losses  [23]. 

£+  and  £_  satisfy  the  boundary  conditions  £+=r£_  at  the 
Z= 0  boundary  and  rE+=E_  at  the  Z=L  boundary  ( L  is  the 
cavity  length  and  r~  0.53  is  the  reflection  coefficient).  The 
other  quantities  in  Eq.  (1)  are  constants:  k ,  A/,  and  T  are  the 
wave  number  (in  the  material)  associated  with  the  resonance 
optical  frequency,  the  electron  density  in  the  active  region, 
and  the  overlap  factor  between  the  laser  mode  and  the  active 
region,  respectively. 

Figure  2  shows  spectra  that  were  obtained  by  solving  nu¬ 
merically  Eqs.  (1)  with  the  following  parameters:  for  the 
saturable  absorber,  we  used  y=  1  0-8  obtained  from  two- 
dimensional  mode  simulations,  assuming  a  n2~  10-9  ^7- 
[23].  The  index  change  due  to  this  n2  at  typical  intracavity 
intensities  is  about  10-3.  The  other  parameters  are  €0 
=  5  cm-1,  T]  =  0.5  ps  [19],  T 2= 0.07  ps  (corresponding  to  a 
gain  full  width  at  half  maximum  bandwidth  of  4.8  THz),  L 
=  0.3  cm,  and  n  =  3,  which  are  typical  values  for  these  lasers. 
N  and  T  are  not  needed  as  long  as  the  pumping  is  expressed 
relative  to  the  lasing  threshold.  Note  that  the  simulated  spec¬ 
tra  presented  in  Fig.  2  are  averaged  over  about  a  microsec¬ 
ond.  Only  then  does  the  average  spectrum  reach  a  steady 
state  and  a  clear  pattern  shows  up.  The  averaging  is  moti¬ 
vated  by  the  fact  that  experimentally  the  spectra  are  acquired 
over  an  even  much  longer  time  scale.  The  envelopes  of  the 
spectra  show  two  clear  peaks  whose  separation  compares 
well  with  twice  the  Rabi  frequency,  similarly  to  the  experi¬ 
ment. 

The  lowering  of  the  RNGH  instability  threshold  by  a  satu¬ 
rable  absorber  can  be  established  analytically  by  means  of 
linear  stability  analysis.  We  propose  this  mechanism  as  the 
main  reason  for  the  observation  of  the  RNGH  instability  at 
much  lower  pumping  than  RNGH  theory  predicts.  In  order  to 
support  this  idea,  we  now  present  spectra  from  another  de¬ 
vice  similar  to  the  one  described  previously.  The  only  differ¬ 
ence  between  the  two  lasers  is  a  shorter  optical  wavelength 
(5.25  71111)  (wafer  no.  3252)  and  a  wider  active  region 
(5  /7,m).  The  two-dimensional  waveguide  simulations  indi¬ 
cate  a  much  weaker  Kerr-lensing  effect  in  these  QCLs  (y  is 
smaller  by  a  factor  of  4),  due  to  the  much  larger  ratio  of 
active  region  width  to  wavelength.  The  measured  optical 
spectra  obtained  at  300  K  in  cw  mode  with  the  5-/7,m  device 
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FIG.  3.  Optical  spectra  vs  pumping  ratio  above  threshold  ob¬ 
tained  in  cw  at  300  K  with  a  5-/um-wide  buried  heterostructure 
lasers  emitting  at  5.25  /um. 

are  shown  in  Fig.  3.  The  data  clearly  show  that  the  laser  is  at 
first  single  mode  close  to  threshold  and  becomes  multimode 
immediately  after  a  slight  increase  of  the  pumping  current. 
The  envelopes  of  the  spectra  consist  of  multiple  peaks,  with 
an  average  separation  0.2  THz,  independently  of  the  pump¬ 
ing.  Numerical  integration  of  Eq.  (1)  without  a  saturable  ab¬ 
sorber  (y=0)  leads  to  spectra  that  qualitatively  agree  with 
the  ones  in  Fig.  3. 

Reference  [15]  suggested  that  the  suppression  of  the  cen¬ 
tral  peak  in  RNGH-type  spectra  can  be  due  to  the  complex 
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level  structure  of  the  gain  medium,  a  dye  molecule  in  that 
case.  We  show  that  SHB  can  also  result  in  the  suppression  of 
the  central  peak  (Fig.  2). 

Our  postulation  of  saturable  absorption  due  to  Kerr- 
lensing  is  supported  by  more  extensive  study  of  different 
devices  beyond  those  shown  in  this  Rapid  Communication. 
First,  we  observed  that  for  the  same  emission  wavelength,  a 
broad  active  region  leads  to  a  less  pronounced  RNGH-type 
signature.  Second,  we  have  also  tested  several  standard  ridge 
waveguide  QCLs,  for  which  the  sidewalls  of  the  ridges  are 
covered  by  the  gold  contact.  For  these  devices  the  coupling 
between  the  optical  mode  and  the  metal  is  expected  to  be 
stronger  and  so  is  the  effect  of  saturable  absorber  due  to  Kerr 
lensing.  The  spectral  behavior  observed  in  this  class  of  de¬ 
vices  is  dominated  by  RNGH-type  instability. 

In  summary,  a  coherent  multimode  instability  in  quantum 
cascade  lasers  (QCLs)  has  been  observed.  It  is  similar 
in  many  ways  to  the  Risken-Nummedal-Graham-Haken 
(RNGH)  instability.  The  threshold  of  the  observed  phenom¬ 
enon  is  significantly  lower  than  in  the  original  RNGH  insta¬ 
bility,  which  is  attributed  to  the  presence  of  a  saturable  ab¬ 
sorber  in  the  laser.  For  devices  with  a  weaker  saturable 
absorber,  the  envelope  of  the  optical  spectrum  consists  of 
many  maxima  whose  separations  are  independent  of  the  int¬ 
racavity  power.  The  nontrivial  shape  of  the  spectrum  can  be 
explained  by  SHB. 
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